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ABSTRACT 

The origin of the diffuse hard X-ray (2-10 keV) emission from starburst galaxies is a long-standing problem. 
We suggest that synchrotron emission of 10 - 100 TeV electrons and positrons (e^) can contribute to this 
emission, because starbursts have strong magnetic fields. We consider three sources of e^ 1 at these energies: 
(1) primary electrons directly accelerated by supernova remnants; (2) pionic secondary e^ created by inelastic 
collisions between CR protons and gas nuclei in the dense ISMs of starbursts; (3) pair e ± produced between the 
interactions between 10 - 100 TeV 7-rays and the intense far-infrared (FIR) radiation fields of starbursts. We 
create one-zone steady-state models of the CR population in the Galactic Center (R < 112 pc), NGC 253, M82, 
and Arp 220's nuclei, assuming a power law injection spectrum for electrons and protons. We consider different 
injection spectral slopes, magnetic field strengths, CR acceleration efficiencies, and diffusive escape times, and 
include advective escape, radiative cooling processes, and secondary and pair e . We compare these models to 
extant radio and GeV and TeV 7-ray data for these starbursts, and calculate the diffuse synchrotron X-ray and 
Inverse Compton (IC) luminosities of these starbursts in the models which satisfy multiwavelength constraints. 
If the primary electron spectrum extends to ~ PeV energies and has a proton/electron injection ratio similar 
to the Galactic value, we find that synchrotron contributes 2 - 20% of their unresolved, diffuse hard X-ray 
emission. However, there is great uncertainty in this conclusion because of the limited information on the CR 
electron spectrum at these high energies. Inverse Compton emission is likewise a minority of the unresolved 
X-ray emission in these starbursts, from 0. 1 % in the Galactic Center to 10% in Arp 220's nuclei, with the main 
uncertainty being the starbursts' magnetic field. We also model generic starbursts, including submillimeter 
galaxies, in the context of the FIR-X-ray relation, finding that anywhere between and 16% of the total hard 
X-ray emission is synchrotron for different parameters, and up to 2% in the densest starbursts assuming a E~ 2 - 2 
injection spectrum and a diffusive escape time of 10 Myr (E/3 GeV)" 1//2 (/i/100 pc). Neutrino observations 
by IceCube and TeV 7-ray data from HESS, VERITAS, and CTA can further constrain the synchrotron X-ray 
emission of starbursts. Our models do not constrain the possibility of hard, second components of primary e 
from sources like pulsars in starbursts, which could enhance the synchrotron X-ray emission further. 

Subject headings: cosmic rays - galaxies: starburst - X-rays: galaxies - gamma rays: galaxies - radio contin- 
uum: galaxies - radiation mechanisms: nonthermal 



1. INTRODUCTION 

Starburst galaxies are intense generators of cosmic rays 
(CRs), which are accelerated by supernova remnants or 
other star-formation processes. CR protons in starbursts 
can produce 7-rays, neutrinos, and secondary electrons and 
positrons. Whatever their origin, CR electrons and positrons 
can produce electromagnetic emission across the spectrum: 
bremsstrahlung losses produce 7-rays; Inverse Compton scat- 
tering of ambient photons produces a broadband spectrum ex- 
tending into 7-rays; synchrotron emission is responsible for 
the non-thermal GHz radio emission. 

Starburst galaxies are observed to be luminous in hard 
X-rays (here defined as ~ 2 — 10 keV) as well. The to- 
tal hard X-ray emission from star-formation typically has 
a luminosity that is 10 -4 times that of the bolometric lu- 
minosity of the starb urst, and is somet i mes used as a star- 
formation indicator (iFranceschini et all 120031: iRanalli et alJ 
120031; iGrimm et al.1 120031; iPersic et alJ |2004 . Most of 

1 Janksy Fellow 

- Institute for Advanced Study, Einstein Drive, Princeton, NI 08540 

3 Department of Astronomy, The Ohio State University, 140 West 18th Av- 
enue, Columbus, OH 43210. 

4 Center for Cosmology & Astro-Particle Physics, The Ohio State Univer- 
sity, Columbus, Ohio 432 1 0, US A 

5 Alfred P. Sloan fellow 



these X-rays are from point sources, but there appear to 
be additional diffuse components. The diffuse hard X- 
ray emission is best studied in the nearby starburst M82 
(IStrickland & Heckmanl 120071) . but a similar diffuse hard 
component is obse rved in NGC 253 ({Strickland et al] 120001 ; 
IWeaver et al.|[2002h . The diffuse hard X -ray emission bears 
some superficial res emblance to the Galactic Ridg e emission 
of the Milky Way (IStrickland & Heckmanl 120071) . although 
that h as recently been resolv ed into stellar sources by Chan- 
dra dRevnivtsev et alJ 120091) . At the other extreme, unre- 
solved hard X-ray emission is also observed in brighter star- 
bursts m^^^^^J^a^nausl^^sdXjSilsxy (LIRG) NGC 
3256 dMoran et alJI 19991; IDra et al.ll2002l) and the prototyp- 
ical Ultraluminous Infrared Galaxy (ULIRG ) Arp 220 (e.g., 
IClements etaT]l2002t iMcDowell et aLll2003l) . Especially in 
these more distant galaxies, it is unclear whether these com- 
pone nts are truly diff use or simply faint, unresolved sources 
(e.g.. lLiraet"aT1l2002l) . 

The diffuse hard X-ray emission typically has a power- 
law continuum spectrum, with the possible addition of 
softer thermal em i ssion c omponents and spectra l lines 
(IPersic & Rephaelil l2002t IStrickland & Heckmanl 120071: 
Lehm eret alJ 12010). Thermal emission from hot plasma 
is one possible source of the hard X -rays: hot gas is 
predicted by superwind theories (e.g., IChevalier & Cleggl 
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119851: IStrickland & Heckmanl 120091) . and thermal emission 
is clearly de tected in soft X-rays (e.g., [Ptak et all 119971 : 
Dahlem et all Il998t [Strickland & Stevensl |2000j), although 
Strickland & Stevens! d2000h argue that the soft X-rays come 
from a cooler phase of gas t han the bulk of the superwin d 
(see also the discussion in Stricklan d & Heckmanl [2009). 
The detection of 6.7 keV iron K lines implies the existence 
of hot plasma tha t could be a source of the hard X-rays 
(Persi c etai] fl998t ICappi et alj [19991: llwasawa et ail 120051 
2009), although it is not clear tha t such emission could 
expla in all of the hard continuum (IStrickland & Heckmanl 
2007). The hard X -ray emission is often attributed to 
unresolved X-ray binarie s, particularly high-mass X-ray 
binar ies (HMXBs) (e.g.. IFabbiano et al.1 119821: IFa bbiano 
19891: iGrifnths&Padovanil Il990t iDavid et al.l 119921: 
Persic & Rephaeli 2002: iGrimm et al.l 120031: iPersic et alj 
20041) . HMXBs in the Milky Way and Magellanic 
Clouds have a power law continuum (photon spectra of 
dN/d E oc E~ r ) with a spectral slope of T « 1.2 (IWhite et alJ 
1983). Many starburst ga laxies have total hard X-ra y emis- 
sion with T « 1-1.5 (e.g jRephaeli et al.ll 199 ill 19951) . which 
makes HMXBs an attractive candidate for the source of 
most of the hard X-ray emission. However, in M82, diffuse 
hard X-ray emission remains even after subtracting point 
sources and extrapolating the luminosity function to faint 
luminosities, and the diffuse e mission is softer (IVs ~ 2-3 ) 
than expected from HMXBs (IStrickland & Heckmar]|2007l) . 
unless the HMXB spectrum cuts off at energies < 10 keV. 

Another possible expla nation for this X-ray emission is In- 
verse Compton emissio n (Hargrave 1974; Schaafetal. 1989; 
Moran & Lehnertl 1997t iMoran et al.ll999tlPersic & Re phaeli 
20031) . CR electrons and positrons (e ) are known to be 
present in starbursts from their synchrotron radio emission, 
and the intense infrared emission of starbursts provides many 
target photons to be upscattered to higher energies. The IC 
spectrum is expected to extend down to the X-rays and even 
lower energies. However, recent estimates generally sug- 
gest that IC emission is too weak by a factor of > 10 in 
M82 and NGC 253 to explain the hard X-ray emi ssion (e.g., 
IWeaver et al.1 120021: IStrickland & Heckmanl 120071) . As with 
HMXBs, the spectral slope of the diffuse X-ray emission in 
some starbursts may be difficult to explain with IC. The CR 
e spectrum around ~ 100 MeV and the resultant IC spec- 
trum at keV energies is expected to be hard (T ps 1.0- 1.5), 
whereas the d iffuse X-ray emission is ofte n softer, as is the 
case for M82 ( IStrickland & Heckmanll2007l) . 

Here, we consider an alternative source of hard X-rays: 
synchrotron emission from CR e with energies above a TeV. 
Synchrotron has previously been c onsidered as a source of 
Galactic diffuse X-ray emissi on dProthero e & Wol fendalel 
Il980t iPorter & Protherod \\Wk . but explaining the diffuse 
X-ray emission from the Galaxy requires CR electrons of 
extreme energies (~ 100 TeV), and stronger magnetic fields 
(B > 20 /iG) than in the diffuse ISM to avoid conflict with 
constraints on Inverse Comp t on emiss i on fro m CASA-MIA 
dAharonian & Atovanl 120001: iBi et alj 120091) . Starbursts 
are expected not only to have large CR populations but 
also stronger magnetic fields than the Milky Way. The 
conclusion of strong magnetic fields in starbursts is sup- 
ported by a number of lines of evidence: (1) minimum 
energy estimates applied to radio detections of starbursts 
imply B ss 50-150 tiG (e.g.. IVolk et al J 1 1 989b iBecld l2005h 
IThompson etai]|2006t IPersic & Rephaelill2010t iBecklboi ll) : 
(2) detailed modeling of CR populations in starbursts with 



fitting of the radio and now 7-ray spectra im ply B « 100 /iG 
in the Galactic Center (ICrocker et all 1201 lbl) B sa 200 p,G 
in the nearby starbursts M82 and NGC 253 (Paglion e et alJ 
1996 1 : iDomingo-Santamarfa & Torres! 120051: IPersic et alj 
20081: Ide Cea del Pozo etai] l2009ah iRephaeli et alJ 120091: 
de Cea del Pozo et al.ll2009bl) . and B>\ mG in the ULIRG 
Arp 220 djorresl |2004); (3) considerations of the linear 
FIR-radio correlation that applies to quiescent star-forming 
galaxies and starbursts, and the strong Inverse Compton, 
bremsstrahlung, and ionization losses exp e cted for CR 
electr o ns in these galaxies (e . g„ IVolkl 119891: [Condon et alJ 
199 It IThompson etai] 120061: iMurphyi 120091: lLacki etai] 
2010); (4) measurements of Zeeman splitting in ULIRGs, 
although these apply t o the denser regions of the starbursts 
(Robis haw et al.l 120081) : (5) constraints on leptonic 7-ray 
emissio n from the Galactic Center, limiting the number of 
CR e ± (ICrocker et al.ll20l0ah . 

Strong magnetic fields imply not only more synchrotron 
power per particle, but that lower energy CR electrons can 
produce synchrotron X-rays. The synchrotron emission of a 
CR e ± peaks near vq = 3j 2 eBsina/(4irm e c) for an electron 
with Lorentz factor 7 with a pitch angle a with respec t to a 
magnetic field of strength B ( Rybicki & Lightman 19791) . For 
an isotropic distribution of electrons, (sin a) = it /A, and this 
translates to a synchrotron emission energy of E c = hvc- 

E c k 1.0£? OTeV fi 2 ookeV, (1) 

where E l0Te v =E e /{\0 TeV) and B 2 oo = 5/(200 /iG) is the 
magnetic field strength in the starburst. We see that CR e^ of 
energy in the range 10-30 TeV will emit X-ray synchrotron 
emission in the 1-10 keV range. 

In this paper, we show that under the most optimistic as- 
sumptions, the synchrotron emission can explain the diffuse 
hard X-ray emission of starbursts, especially Arp 220. The 
synchrotron spectrum (vL v ) rises at 1 GHz frequency and 
peaks near 10 - 100 GHz, but remains constant (or slowly 
falling) all the way to X-ray frequencies. There are three main 
reasons for this. (1) The synchrotron losses of GHz-emitting 
e in starbursts are ~ 1 /5 - 1 /20 of t he total losses including 
bremsstrahlung, ionization, and IC (IThompson et al.l 120061 : 
lMurphy||2009l: lLacki et al.ll20 1011201 ll) . Thus little of the en- 
ergy in these CR e is transformed into synchrotron emis- 
sion. However, while ionization and bremsstrahlung losses in- 
crease slowly with CR e^ energy, synchrotron (and IC) losses 
become faster increasing with CR e ± energy. Therefore, 
bremsstrahlung and ionization become c omplet ely unimpor- 
tant at 10 TeV (see the loss time scales in § 12.1.11 ). (2) Further- 
more, IC cooling will be suppressed at energies above 10 TeV. 
Most of the starlight generated by starbursts is efficien tly ab- 
sorbed by dust and re-emitted in the IR (c.f. Volk 1989]), and 
the majority of this IR is blackbody emission from ~ 40-50 K 
dust grains. For FIR photons of energy 6ir and wavelength 
A w 100 /J,m, 20 TeV is near the characteristic energy scale 
for Klein-Nishina suppression of IC cooling: 

2 4 
m c 

£KN«-^«21TeV(A/100^m). (2) 
£ir 

(3) A more subtle but important effect of the ~ 20 TeV Klein- 
Nishina cutoff is that the threshold for pair production for 7- 
rays on the FIR emission is near this energy. Thus the in- 
tense IR emission of starbursts converts the 7-ray emission 
at 10 - 100 TeV energies into 10 - 100 TeV e ± . Since the 
power injected in pionic 7-rays is greater than the power in 
pionic secondary e , these tertiary pair e can dominate the 
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CR spectrum at 10 - 100 TeV if the pair production optical 
depth is greater than unity (77 attenuation and pair production 
at 10 - 100 TeV is expected to be small in galaxies like the 
Milky Way, see | Ma stichiadis et al J 1 1 99 It iMoskalenk o et al.l 
l2006t IStawarz et al.l 120101) . providing an additional popula- 
tion that can radiate X-ray synchrotron. 

We begin in §|2]by presenting order of magnitude estimates 
showing that synchrotron may be important. Then, in later 
sections we construct one-zone models to evaluate these ef- 
fects using standard assumptions for CR modelling and in the 
context of extant multiwavelength data. In § [3] we describe 
our one-zone models of the CR spectra, including the pair 
produced e . In § @] we present the results of these mod- 
els for the Galactic Center, NGC 253, M82, and the nuclear 
starbursts of Arp 220. We discuss some broader implications 
in §[5] including the sync hrot ron contribution to X-rays from 
submillimeter galaxies (§ 15. U and further constraints on syn- 
chrotron X-rays from future neutrino and TeV experiments 
(§ 15.2b . Finally, we conclude by discussing future work that 
can be done in § [6] 

2. MOTIVATION 
2.1. Relevant Cooling Processes 
2.1.1. e ± losses 

Starbursts contain dense gas and strong radiation fields, 
which can cool CR e ± through bremsstrahlung, ionization, 
and IC emission. These cooling processes have associated en- 
ergy loss rates for individual CR e ± &brems> &ion, bic, respec- 
tively, which can be comparable to the synchrotron energy 
loss rate b sync h- In addition, CR e ± can escape by advection 
or diffusion. These processes compete with synchrotron emis- 
sion in the magnetic fields of starbursts for the kinetic energy 
ofaCRe ± . 

The synchrotron cooling time of CR e ± is E /b syriC h(E) ~ 
6wm 2 c 4 /(ca T B 2 E) dRybicki & Lightmanll 19791) . or 



fsynchC^) » 31 yr 



10 TeV 



B 



200 fiG 



(3) 



For energies greater than 1 TeV, and magnetic fields of at 
least 200 /iG, this is about the light crossing time for 100 pc, 
roughly the scale height for a nuclear starburst. Therefore, at 
these energies, e^ are in the "cal orimeter" lim it, in which the 
CR e ± cool before escaping (c.f., Volk 1989). 

For a gas density n, the ionization cooling time of CR e ± is 
E/b ion (E), or 



t lon (E)~ 1.9 x 10 9 yr 



(4) 



10 TeV J V 250 cm" 3 

from e.g., Schlickeiser (2002), where we set the ln7 term to 
its value at 10 TeV. The bremsstrahlung cooling time is 



1.2xl0 5 yr( " _, )"', (5) 
V 250 cm - 1 J 



^bremsO^) ' 

at high energies ( assum ing «e e ~ 0.1/i#; see 
IStrong & Moskalenkol 119981) . Synchrotron losses 

therefore dominate bremsstrahlung losses as long as 
B > 10 nG (E/TeVT 1/2 (n/250 cm" 3 ) 1 / 2 , and syn- 
chrotron losses dominates ionization losses so long 
as B > 0.2 /iG (£/TeV)" I (n/250 cm" 3 ) 1 / 2 . However, 
bremsstrahlung and ionization can dominate at the ~ GeV 
energies where e are responsible for GHz em ission 
( Tho mpson etaf]|2006t lMurphy||2009t lLacki et alj|20!ob . 



The Inverse Compton cooling time is more complex be- 
cause Klein-Nishina effects appear at the relevant ~ 10 — 
100 TeV e nergies in the FIR-dom i nated radiation fields of 
starbursts. ISchlickeiser & Ruppell d2010l) show that the IC 
loss time in a greybody radiation field with temperature T and 
radiation energy density t/ raa can be approximated as 



tic' 



3m e c 2 7I + 7 2 
4ca T U md 77! 



(6) 



where j K m 0.21m e c 2 /(k B T) = 4.0 x 10 7 (r/40 K) _1 . Thus, 
above ~ 20 (T/40 K)" 1 TeV, Klein-Nishina effects will cause 
the IC energy loss time to grow as 7, while the synchrotron 
loss time continues to fall as 7" 1 . If Ub < U m d, then syn- 
chrotron dominates IC when 7 > ^k^V^IUb- 1, or: 



E > 20 TeV 



40 K 



Urud-Ut 



1/2 



(7) 



The nearby, prototypical starbursts in M82 an d NGC 253 
have t otal infrared luminos ities of ~ 10 10 5 L© (Me lo et all 
2002t ISanders et alj l2003h and rad i i of ~ 200 pc (e.g. 
Turner & Holll983l:lGoetz et al|ll990l : IUlvestad & Antonuccil 
19971; IWilhams & Bowerfl2010h . Plugging in these specific 
values for a disk geometry, the radiation field has an energy 
density of approximately U « L/(2irR 2 c), or 



IU~ HOOeV cm" 



R 



1O 1O - 5 L 0/ / V200pc 
This gives us an IC loss time of 



f Thomson , 
'IC 



300 yr 



1O 1O5 L 
in the Thomson regime, and 
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'ic 



:71yr 



L 



R 



, (—)' 
200 pcy v TeV / 



( T 



(8) 



(9) 



200 pc/ V40K 



100 TeV 
(10) 

in the extreme Klein-Nishina limit. For the specific values 
we have been using - L = 10 10 5 L Q , R = 200 pc, T = 40 K - 
equation [7] shows that synchrotron dominates over IC above 
~ 5 TeV when B = 200 ^G and - 35 TeV when B = 100 /xG. 

Therefore unless the magnetic field energy density is much 
lower than the radiation field energy density, synchrotron 
losses will be the dominant loss process at a few tens of TeV. 

2.1.2. Proton losses 

CR protons can modify the e^ population by creating pio- 
nic secondary electrons and positrons through collisions with 
interstellar gas atoms. The time for CR protons to lose all of 
their energy through this process is given as 



*pion = 2x 10 5 yr 



( 1 V 1 

V 250 cm" 3 / 



(11) 



bv lMannheim & Schlickeiserl d!994l) . Different sources in the 
literature give different pionic loss times that can be larger by 
a factor ^2-3; however, these decrease slightly with proton 
energy as the pionic cross section increases (see the discussion 
in Appendix lAl. so that at TeV energies the pionic loss time 
is closer to equation [TT] 

Since protons have much longer cooling times than TeV 
CR e ± , it is not as obvious whether they can escape. Advec- 
tion is clearly present in starbursts in the form of the observed 
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large-scale winds ( e.g. jHeckman et al.ll2000tlHeckmanl2003l; 
IHeesen et al.ll2.Ql II) . The wind crossing time for a starburst of 
scale height h is h/v, or: 



fwind = 3 x 10 5 yr 



( h 



VlOOpc / \300kms- 1 



(12) 



In starbursts such as M82 and NGC 253, the advective and 
pionic lifetimes are expected to be roughly equal. The 
nuclear starbursts of Arp 220, with («) ps 10 4 cm" 3 (e.g., 
iDownes & Solo mon 1998|), have very short pionic loss times 
(~ 5000 yr), indicating that they are "proton calorimeters": 
most of the power injected into CR protons with energies 
above the pion-production threshold is lost through pionic in- 
teractions. 

However, a key uncertainty is whether CRs sample gas of 
the average density. The 7-ray luminosity of the Galactic 
Center region (R < 112 pc) relative to its star-formation rate 
indicates this is not the case for that region, as does the ra- 
tio of synchrotron ra di o from CR e to inf rared emission 
dCrocker et aU201 1 alibi) . ICrocker et ail d2011bl) explains these 
observations as being caused by a powerful wind in the Galac- 
tic Center region advecting CRs out of the disk before they 
can enter the molecular clouds containing most of the gas. 
Thus CRs experience gas of much lower density than aver- 
age. While the Galactic Center is underluminous in 7-rays 
and radio, this is not true for the starbursts M82 and NGC 
253, which fall on the FIR-radio correlation and have a larger 
7-ray to star-formation ratio than the Milky Way, consistent 
with CRs experiencing average gas densities in these star- 
bursts dLacki et alj|201 lh . 

A final uncertainty is whether diffusive escape plays any 
role. Diffusive escape in the Milky Way has an en- 
ergy dependence t&g oc E~° 3 - E~° 6 that steepens the CR 
proton spectrum, since it is the dominant timescale (e.g., 
iGinzburg & Ptuskir] 119761) . In contrast, the relatively flat 
GeV-TeV spectra of observed starburst regions (M82, NGC 
253, and the Galactic Center) indicate that up to TeV ener- 
gies, an energy-independent process must determine the life- 
times of CR protons; however, at the still higher energies we 
are considering, diffusion can dominate. We consider several 
values of the diffusive escape time to address this uncertainty. 

2.2. Primary Electrons 

Primary electrons dominate the GHz-emitting electron 
population in normal galaxies. However, they are gen- 
erally expected to be sub-dominant in starbursts with re- 
spect t o pionic secondary e^, though still a significant mi- 
nority dRengaraianl 120051: iThompson. Ouataert. & Waxmanl 
120071 lLacki et alj|2010l) . Detailed modelling of CR popu- 
lations in starburst galaxies find that primary electrons are 
subdominant at GeV energies in M82 (|de Ce a del Poz o et all 
' 2009ah. NGC 253 dDomingo-Santamarfa & Torres! 1200$ 
Repha eliet all 120091). an d Arp 220 dT orresll2004l) . although 



models by Crocker e t alJ (1201 lbl) find primary electrons dom- 
inate in the Galactic Center. At higher energies, where the 
protons experience stronger diffusive losses, these models in- 
dicate higher primary fractions. 

It is not clear where the primary electron injection spectrum 
ends in starbursts. As the electrons are accelerated to higher 
energies, they also experience more severe synchrotron and 
Inverse Compton cooling. At some point, the cooling losses 
balance the rate of acceleration, and there can be no further 
acceleration of primary electrons. In the standard supernova 



acceleration theory, equilibrium between cooling and acceler- 
ation occurs when: 



i 27 TeV VscaaB. 



-1/2 



5000« SNR ,200> 



(13) 



where 5OOOV5000 km s" 1 is the speed of the supernova shock 
and 20Qj?sNR.2 00 is the supernova remnant magnetic field 
(Gaisser 1990). From equation[U the end of the primary syn- 
chrotron spectrum will then be at: 



rsynch ^71 ^ y 2 

£ my ~ / . J JVC V V* 



5000 
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\BsnrJ 



(14) 



if the magnetic field strength in the starburst B is similar to 
the magnetic field strength in the supernova remnant. X- 
ray observations of supernova remnants in the Milky Way 
have revealed synchrotron emission from 10 - 100 TeV 
electrons, confirmin g these energies are reached in super- 
nova remnants (e.g., iKovama et alj|1995t lAllen et all 1 1997b 
[Reynolds & Keohand [19991: iRevnoldsl 120081: iRevnolds et all 
2008). Thus, SNRs might be able to accelerate primary elec- 
trons to the energies where they will produce hard synchrotron 
X-rays in the diffuse ISM of starburst galaxies. 

However, it is not clear that SNRs are responsible for all 
of the primary CR electrons; other objects such as superbub- 
bles or pulsars may contribute (e.g., Butt 2009). Pulsar Wind 
Nebulae (PWNe) powered by pulsar spindown after their birth 
supernova in particular may inject a very hard component of 
e ± domin ating at TeV energies, possibly exte nding to 10 - 
100 TeV dYuksel et al.l2009tlBamba et alJ2010h . Primary CR 
e ± from PWNe have been invoked to explain anomalies in 
the CR electron spectrum observed at Earth. The birthrate of 
pulsars is correlated with star-formation, and pulsars should 
be present in large numbers in starbursts (c.f. iPerna & Stella! 
■20041 Mannhei m et alJiMoh . Their spin-down luminosities 
could provide enou gh power to be c omparable to the main 
component of CRs (Lacki et al. 20TTT). 

If a hard primary CR e ± spectrum extends to 10 - 100 
TeV energies or beyond, then the synchrotron X-ray luminos- 
ity may be very bright. Suppose that for every supernova, 
roug hly £ e CR = £48 10 48 erg of CR electr ons are accelerated. 
This could occur if 0. 1 % of a supernova remnant's kinetic en- 
ergy (10 51 ergs) were converted into primary CR electrons, 
which is about 1% of the total energy expected to go into CRs 
dStrong et al. 2010, the exact ratio depends on the poorly con- 
strained low energy CR e spectrum). For an E~ 2 injection 
spectrum extending from 1 GeV to 1 PeV, then E 2 dQ/dE = 
££ R /ln(10 ). At very high energies, a fraction / sync h will go 
into synchrotron, where / synch = t hfe /t synch ~ 1/(1 +4ynchAic): 



/synch 



1 + 



t/rad/t/B 



(15) 



1 + (7/7k) 2 _ 

(eqs. [3] and [6}, which we have argued to be near 1. Since 
the characteristic synchrotron emission energy Eq is propor- 
tional to \flEe), the synchrotron emission of the CR e pop- 
ulation in one dex of electron energy will be spread over two 
dex in synchrotron fre quency. We finally hav e vL v (keN) = 
{\/2)f mch EldQ e /dE e dLoeb & Waxmanll2006l) . or: 



i/L„(keV) ■ 



1.1 X 10 ^48/synch ( —[ ) ergS S 



r 



SN 



(16) 



where r$ N is the supernova rate. For example, M82 is 
believed to have a supernova rate of roughly ^0.1 yr" 1 
(although with large uncertainties), suggesting that its syn- 
chrotron X-ray luminosity from pionic e ± may be 1 x 



10 38 ergs s about 3% of the observed luminosity of 
the diffuse X-ray emission (L (2-8 keV) = (4.4 ± 0.2) x 
10 39 ergs s" 1 , or vL v « 3.2 x 10 39 ergs s" 1 per bin in In en- 
ergy; IStrickland & Heckmarj|2007l) . Softer CR e ± injection 
spectra will have still lower synchrotron X-ray luminosities. 

If SNRs or other discrete sources do accelerate ~ 10 TeV 
primary electrons, the synchrotron X-ray emission may not 
be spread continuously throughout the starburst, but concen- 
trated near the CR sources, because the synchrotron cooling 
time for these electrons is so short (eq. [3}. Even if such elec- 
trons free-stream, they will not travel farther than ~ 10 pc 
from their sources. If there are a small number of CR e^ ac- 
celerators, then the synchrotron X-ray emission should come 
from a few small diffuse regions, just as patchy "c ells" of high 
energy are predicted for the Milky W ay (c.f., Shenl ll970l: 
lAharonian et al.|[l995HAtoyan et al.lll995l) . If starbursts, with 
a high star-formation rate concentrated into a small volume, 
instead contain many accelerators, the primary e^ confine- 
ment regions around these accelerators will overlap and the 
X-ray emission will arise throughout the starburst. Unlike pi- 
onic secondary e or pair production e ± from pionic 7-rays, 
both of which depend on the CR proton spectrum, primary e 
at these energies are unaffected by the escape because cooling 
is so quick. 

2.3. Pionic Secondaries 

CR protons can inelastically scatter off protons in the ISM 
to produce pions, which decay into secondary e^, 7-rays, and 
neutrinos. From the lifetime calculations given in § 12.1.21 
starbursts are expected to convert much more of their CR pro- 
ton energy into pionic products than the Milky Way and ap- 
proach the "proton calorimeter" li mit (Loeb & Waxman 2006; 
Thom pson. Ouataert. & Waxmanl l2007b lLacki etal.1 2010). 
From the HESS detection. lAcero et al.1 (120091) claimed that 
5% of the proton energy was converted into pionic prod- 
ucts in NGC 253, although they assumed a hard 7-ray spec- 
trum. From tins Fermi, HESS, and VERITAS detections, 
lLacki et ail (1201 ll) inferred a proton calorimetry fraction of 
about 1/3 for NGC 253 and M82 from the ratio of the > 
GeV 7-ray and IR luminosities. The secondary e^ are ex- 
pected to dominate in starbu rsts at GeV energies based on 
physi c al considerations (|Reng araian 20 051 lLoeb & Waxman 
2006; [Thompson. Ouataert. & Waxman ~ l2007t lLacki et al. 
20101) and detailed models of M82, NGC 253, and 



Arp 220 (lTorresll2004l: IDomingo-Santamaria & Torresl 120051; 
Ide Cea del Pozo et al.l2009al;lRephaeli et al.l2009t) . Primaries 
probably dominate in the Galactic Center due to a fast wind 
advecting CR protons from th e region before they can interact 
with the molecular cloud gas (ICrocker et al.ll201 lbt) . 

In the Milky Way, primary CR protons (and nuclei) are ac- 
celerated with a power law spectrum extending to energies 
of at least ~ PeV, the so-called "knee" in the CR spectrum. 
Above these energies, the CR spectrum steepens: this may be 
because Galactic CRs are not accelerated to higher energies 
and we are seeing the transition to a different component of 
CRs, or because Galactic CRs propagate differently above the 
knee. Essentially nothing is known about the knee in star- 
bursts, but presumably it is also at least at a PeV in energy. 

If most of the energy injected into CR protons is lost 
to pion production for CR protons of up to a PeV in en- 
ergy, then the large population of secondary pionic 10 - 100 
TeV e produced by pion decay can emit bright hard X- 
ray emission. Suppose that roughly £cr = 10 erg of CR 
protons are accelerated per supernova (that is, roughly 1/10 



of the SN kinetic power goes into CRs; Strong et al. 2010). 
Once again we have E 2 dQ/dE = £cr/1h(10 6 ) for an E~ 2 in- 
jection spectrum extending from 1 GeV to 1 PeV. A frac- 
tion F C!l \(E p ) of that power is lost to pions for a CR pro- 
ton energy E p ; of that, ~ 1/6 will go into secondary e^ 
(e.g., ISteigman & Strittmatterll 19711: lLoeb & Waxmanll2006l) 6 
and / S y nc h of the secondary e^ 1 power in turn goes into syn- 
chrotron emission. From a similar argument as the pri- 
maries, we finally have i/L^ jkeV) = (1 /2)f syn c hE 2 dQ e /dE e = 
(F^/l2)f synch E 2 ,dQ p /dE p (lLoeb & Waxm anl l2006l) . or: 

^(keV) w 1.9 x 10 4() F cal (100 TeV)/ synch f ergs s"\ 

(17) 

where Tsn is the supernova rate. For a supernova rate in 
M82 of 0.1 yr" 1 , we get 2 x 10 39 ergs s" 1 for / synch = F ca i = 1, 
more tha n half the observed luminosity of the diffuse X-ray 
emission (IStrickland & Heckm an 2007). In practice, energy- 
dependent diffusive escape of the primary protons, advective 
escape in the starburst superwind, and softer CR proton injec- 
tion spectra will reduce the secondary e ± luminosity. 

In contrast with primary electrons, pionic e ± will be gen- 
erated throughout the starburst, instead of being concentrated 
near the source of CR protons, because CR protons can travel 
a longer distance during their lifetime. Therefore the syn- 
chrotron X-ray emission secondary e ± will not have as patchy 
a structure as the primaries, but they will be generated wher- 
ever there is gas being sampled by CRs. 

2.4. Pair Production Tertiaries 

Starburst galaxies are generally predicted to become 
opaque to 7-rays above a few TeV, because 7- rays will pair 
produce with the IR light in the starbursts dTorresi l2004t 
IDomingo-Santamaria & Torresll2005l ; Ilnouell201 ll) . Thus, the 
pionic 7-ray emission at 10 - 100 TeV will be efficiently con- 
verted into 10 - 100 TeV e ± . What is less appreciated is that 
these e ± will lose much of their energy to synchrotron cool- 
ing, emitting mostly X-rays. Thus, if multi-TeV CR protons 
in starbursts have efficient pionic losses, then up to half of the 
CR proton energy at these energies will go into synchrotron 
X-rays, with a similar fraction escaping as neutrinos. 

The optical depth of a starburst to multi-TeV photons is 
t 77 rj tent<7 77 , where h is the height of the starburst disk, cr 77 
is the cross section of pair-production, and «ir is the number 
density of target IR photons. Near threshold (where E 1 s» £kn 
from eq. [2]», er 77 sa 0-7-/4. The average energy of IR photons 
in a greybody field is 7r 4 £(3)AT/30 ~ 2.10kT, where £ is the 
Riemann zeta function. The number density of IR photons 
for a starburst disk of radius R (and emitting area 2irR 2 ) is 
approximately «ir rj L\^/{2hR 2 c x UOksT), where T is the 
typical temperature of the IR photon. We have 



'77 



5.8 



-FIR 



R 



100 pc J V200pc J V 40K 



/ T 



(18) 

demonstrating that luminous starbursts are opaque to ~ 
30 TeV photons, turning them into pair e . In practice, since 
about half of the total infrared radiation is FIR (Calzet tTet al.1 

6 This is because the charged pions receive 2/3 of the energy, and each 
charged pion ultimately decays into a positron or electron and three neutrinos 
of roughly equal energy. Therefore, recieve 1/4x2/3 = 1/6 of the pion 
energy. Similarly, neutrinos receive 1/2 of the pion energy, and 7-rays recieve 
the remaining 1/3, from the neutral pions. 
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2000), the expected r 77 will be lower by a factor of a few, 
but still of order unity. This is in contrast t o the Milky Way, 
which is essentially transparent to 7-ray s (Mastichiad is et alj 
U99UlMoskalenko et al.ll2006tlStawarz et al.ll2010l) . 

Pionic 7-rays are expected to dominate the VHE 7-ray lu- 
minosity of starbursts, so the calculation of the synchrotron 
power from tertiary pair e ± is similar to that for secondary 
pionic e . CR protons roughly inject 2 times more energy in 
pionic 7-rays than in pionic secondary e^, so the synchrotron 
power should likewise be twice as great for pair e ± than direct 
secondaries (eq.fTTIi: 



i/L„(keV) « 3.8 x 10 40 F ca i(100 TeV)/ s - 



ynch 



yr 



SN \ -1 

— ) ergs s 



(19) 

for an £ 2 injection spectrum and 10 50 ergs per supernova 
in CR protons. Again comparing to M82 with a supernova 
rate of Tsn ~ 0.1 yr , we find that vL v may be as high as 
4 x 10 39 ergs s" 1 (with F ca i = / sy nch =1), equal to the observed 
diffuse X-ray emission. As with secondary pionic e ± , diffu- 
sive and advective losses (lower F ca \) and softer injection spec- 
tra reduce the predicted synchrotron X-ray luminosity. On the 
other hand, non-pionic 7-rays, such as from discrete sources 
like pulsars, will also produce pairs in the starburst radiation 
field, and this can enhance the electron population and syn- 
chrotron X-ray luminosity further. 

Like secondary e^, tertiary will not be concentrated 
near CR accelerators, since the pionic 7-rays that generate 
them are emitted everywhere in the starburst region. How- 
ever, the pair e ± production may be enhanced near luminous 
IR sources within the starbursts. A treatment of this effect 
requires a radiative transfer calculation, which is beyond the 
scope of this work. 

3. MODELLING ASSUMPTIONS 

To understand the synchrotron and Inverse Compton X- 
ray emission of starburst galaxies, we create models of the 
steady-state CR populations of the Galactic Center, NGC 
253, M82, and the nuclei of Arp 220. Our goal here is to 
sketch out the parameter space allowed by multiwavelength 
data using a few standard assumptions, and investigate the 
synchrotron and IC emission that arises under these assump- 
tions. We model the starbursts as one-zone disks of radius R 
and midplane-to-edge scale heights h. The evolution of CR 
population is governed by the diffusion-los s equation (e.g., 
iGinzburg & Ptuskinlll976l:IStrong et al.ll2007l) . In steady-state 
one-zone models with no spatial or temporal dependence, the 
diffusion-loss equation reduces to the leaky box equation: 

N(E) d 

■—[b(E)N(E)]-Q(E)=0. 



... /r (20) 

t Me (E) dE 

Here, Q(E) is the injection spectrum of CRs, b(E) is the cool- 
ing rate of CRs (including ionization, bremsstrahlung, IC, 
and synchrotron), and ?Hf e (£) is the lifetime of CRs from es- 
cape (both diffusive and advective) and pionic losses for pro- 
tons. The CR lifetime fufeCE) is calculated as [tm e (E)] = 
[tfoffiE)]' 1 +t~fa + [t 7r (E)]~ 1 , where fdiff(E) is the diffusive es- 
cape time, / a dv is the advective escape time, and t^{E) is the 
pionic los s time for protons . We use the numerical code de- 
scribed in lLacki et al.l (1201 Oh to find the st eady-state CR spec- 
tra, employing a Green's function given in Itorresl (120041) . See 
Table [T]for the parameters we used for each starburst. 

3.1. Injection 



We calculate the star-formation rate by directly converting 
the TIR (total infrared; 8- 1000 /im) luminosity of the star- 
burst disk: 

SFR = L TIR /( £C 2 ), (21) 

where £ = 3.8 x 10 -4 is a dimensionless factor relating the lu- 
minosity to the instantaneo us star formation rate, and is IMF 
dependent (Kennicutt 1998). The factor e is derived assuming 
a starb urst that is cont inuous over 10 - 100 Myr and a Salpeter 
IMF (lKennicutjll998h . The CR energy injection rate per unit 
volume is assumed to be proportional to star-formation rate: 

e CR .e = 9.2x 10- 6 £5i^i7«£/0.001)/V (22) 
e CR ,p = 9.2x 10- 4 £ 5 iV'i7W'7/0.1)/V (23) 
= e CR ,e<5 (24) 

(25) 

for electrons and protons respectively, where £ is the electron 
acceleration efficiency, S is the ratio of proton accleration effi- 
ciency 77 to £, i/in = (/?sn/£)/(17 Mq ), V = 2nR 2 h is the star- 
burst volume, and /?sn is the SN rate per unit star formation. 
The total SN rate in the starbursts is r SN = 17L T ir/(M c 2 ): 



r SN = 0.036 yrVn 



10 



10.5 



(26) 



We assume ^17 = 1 throughout this work. The efficiency of 
primary CR electron acceleration is described by the £ pa- 
rameter, and the ratio of energy going into CR protons and 
CR electrons is 8. 

We assume that CR protons and electrons are respectively 
injected with a momentum power law spectrum 

dQ p /dq = C p q- p (27) 
dQ e /dq = C e q- p (28) 

per unit volume, where q is the CR proton or electron momen- 
tum. A minimum kinetic energy cutoff of K m \ n = 1 MeV was 
used with these spectra. The CR proton spectrum is assumed 
to extend to 7 max = 10 6 , corresponding to an energy of 938 
TeV (7max = 10 5 and 10 7 are considered in Appendix IB! . We 
try different cutoffs in the Lorentz factor 7m™ of the primary 
electron spectrum, both 10 6 and 10 9 . We note that we are as- 
suming the "test particle" approach to CR acceleration here. 
Nonlinear effects can result in more complicated spectra, with 
breaks in the power l aw at low energy and spec t ral hardenin 
at high energy (e.g.. iBerezhko & Ellison! [19 99: Ellis on et 
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2000; Mal kov & O'C Druryl200lHBlasi et alJl2005l) . We also 
neglect additional, harder components of from different 
CR accelerators. In general, a hardening of the injected elec- 
tron spectrum at high energy will result in more synchrotron 
X-ray emission and a softening of the injected electron spec- 
trum will result in less synchrotron X-ray emission. 

The normalizations of the CR injection spectra are then set 
by calculating the integral of the kinetic energy injected and 
equating with ecR- 



9max.p 



£cr,p = C p I 1 P [\ < ? 2c2 + m P c4 ~ m p c ) dq (29) 
£CR.e = C e / q~ p {^q 1 c 1 + rn 2 e c A -m e c 7 \ dq (30) 



where q min = (l/c)JK^ in + 2K min mc 2 and q mdx = 



racy (7m™) 2 - 1. The ratio of energy in CR electrons to 
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TABLE 1 
Model Parameters 



Parameter Units Galactic Center NGC 253 Core M82 Arp 220 West Arp 220 East 



Assumed parameters 



B 


fiG 


50 - 100 


50 - 400 


50 - 400 


250 - 16000 


250 - 16000 




g cm -2 


0.003-0.1 


0.10 


0.17 


10 


10 


h 


pc 


42 


50 


100 


50 


50 


R 


pc 


112 


150 


250 


50 


50 




pc 


112 


150 


250 


50 


50 


D 


Mpc 


0.008 


3.5 


3.6 


79.9 


79.9 


^ ? wind 


km s 


600 


300 


300 


300 


300 


^TIR 


L.Q 


4 x 10 8 


2x 10 10 


5.9 x 10 10 


3 x 10" 


3 x 10 11 


Lx (diffuse) b 


ergs s 1 


7.4 X 10 36c 


8.5 X 10 38d 


4.4 x 10 39c 


4 x 10 40f 


1.5 x 10 40f 



Fiducial Parameters 



B 


/iG 


100 


100 


150 


4000 


4000 


s s 


g cirT 2 


0.003 


0.10 


0.17 


10 


10 


P 




2.2 


2.2 


2.2 


2.2 


2.2 


tm(3 GeV) 


Myr 


1 


1 


10 


10 


10 


prim 

/max 




to 9 


10 9 


10 9 


io 9 


10" 


V 




0.14 


0.40 


0.10 


0.10 


0.10 






0.0066 


0.027 


0.0077 


0.021 


0.0083 


Smermd.OGHz) 


Jy 


370 


0.63 


1.3 







Fiducial Results 



Lx (synch) g ergs s 1 


4.3 x 10 35 


7.2 X 10 37 


9.6 x 10 37 


6.0 x IO 39 


5.2 X IO 39 


r (synch)S 


2.12 


2.08 


2.15 


2.12 


2.12 


L x (IC)S ergs s" 1 


9.0 x 10 33 


2.7 X 10 37 


4.1 x 10 37 


2.6 x IO 39 


1.9 x 10 39 


T (IC)S 


1.58 


1.35 


1.31 


1.44 


1.44 


5 


84 


59 


51 


18 


48 




0.0037 


0.15 


0.33 


0.97 


0.97 


F ca l(> lOTeV) ■■■ 


6 x ur 4 


0.016 


0.11 


0.93 


0.93 



d Adopted sightline distance through starburst to Earth, for computing observed TeV 7-ray spectrum. 

Adopted diffuse hard X-ray emission, with which synchrotron and IC are compared. 
c Galactic Center luminosity in 2 - 10 keV X-ray band from Koyama et alj 119961) , as extrapolated to \£\ < 0.8° and \b\ < 0.3° 
assuming a constant surface brightness. 

^ NGC 253 disk diffuse 2 - 10 keV X-ray luminosity from Bauer et al. 1 2008). Note that this includes the outlying regions of the 
galaxy and not just the starburst core. 

e Luminosity of diffuse hard X-ray excess in 2 - 8 keV band from Strickland & Heckman 1 2007). The uncertainty is 0.2 X 
10 39 ergs S -'. 

Absorption-corrected Arp 220 X-ray luminosities from Clements et al. 1 2002). We assume that Arp 220 X-l is the western 
nucleus and Arp 220 X-4 is the eastern nucleus. 

g Synchrotron and IC luminosities and photon indexes, for the hard X-ray energy bands given for Lx (diffuse) (2-8 keV for M82, 
2-10 keV for the other starbursts). 



CR protons injected at high energies can be approximated as: 

2 \ p~ 2 

S. (31) 



'-I 



In the Milky Way, S « 50-100, which is expected from charge 
conservation in the acceleration reg ion for apw 2.2 m omen- 
tum power law injection spectrum ( Schlickeiser 2002). This 
value of proton/electron injection ratio is also inferred from 
propagation studies that compare with o bservations of CRs 
in the Milky Way (e.g., IStrong et al .1120101) , and from the CR 
pressure derived from shock structure in Tycho's supernova 
remnant dWarren et al.ll2005l) . It is also the approximate ratio 
of the CR proton and elec tron energy densities at Earth (e.g., 
iGinzburg & Ptuskinll 19761) . 

3.2. Propagation 

In solving the leaky box equation for CR e ± , we con- 
sider energy losses from ionization, bremsstrahlung, syn- 
chrotron, and IC. CR protons experience continuous ioniza- 
tion losses, which are important at low energies where pionic 
losses are negligible. All CRs can also escape, either through 
energy-dependent diffusion, or energy-independent advection 



(winds). We assume an E~ x l 2 energy dependence to the diffu- 
sive escape times; the exact energy dependence is not known, 
although our assumption is conservative in that CR protons 
will tend to easily escape at higher energy without producing 
pionic e ± and 7-rays. The normalization of the diffusive es- 
cape time in starbursts is also not known, so we consider a 
variety of normalizations: 



fdiff = fdiff(3 GeV) 



3 GeV 



-1/2 



(32) 



with fdiff (3 GeV) ranging from 1 Myr to 00 (in the Milky Way , 
f di ff(3 GeV) re 30 Myr; IConnelll fl998l IWebber et al.l[2003h . 
Different energy dependences would alter the CR proton pop- 
ulation at TeV energies, and we naturally expect more sec- 
ondary and pair e ± (and their accompanying synchrotron X- 
rays) if there is a weaker energy dependence in fdiff than in 
equation [32] However, the TeV 7-ray luminosity will also 
be different, so in starbursts with a TeV 7-ray luminosity 
constraint, the allowed parameter space would change. This 
would weaken the effect of different energy dependence in 
fdiff, since there is a shorter "lever arm" between the popu- 
lation traced by TeV 7-rays (~ 10-100 TeV protons for 1 - 
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10 TeV 7-rays) and synchrotron X-rays (~ 100- 1000 TeV 
protons). For E~ - 3 energy dependence (the weakest usually 
expected), the diffusion time is greater than 300 kyr at PeV 
energies as long as f<jjff(3 GeV) > 10 Myr, meaning that ad- 
vection (and any pionic losses) will be more important than 
diffusive escape (see equation fT2b. Thus, relatively slow dif- 
fusive escape with a E~ 03 energy dependence is equivalent 
to no diffusive escape at all (fdiff(3 GeV) — > oo). In addition, 
in models where primary electrons dominate the synchrotron 
X-ray flux (the Galactic Center and low B models of other 
starbursts), diffusive escape will have little effect, since the 
synchrotron cooling time is so short at these energies (equa- 
tion rj. 

Advective escape times are more well known; a wind with 
a speed of a few hundred kilometers per second will carry 
a CR out of the starburst in a few hundred kyr (eq [TZb . 
Our models of NGC 253, M82, and Arp 220 assume that 
Vwind = 300 km s" 1 . We use a wind speed of 6 00 km s" 1 for 
the Galactic Center models, as ICrocker et al.1 (1201 lbl) infers 
for the Galactic Center region. 

To calculate the effects of pionic losses on the CR proton 
spectrum, we directly integrate all of the energy going into 
pionic secondary products: 



-K„ 



EfE„ 
/ E s 
+ Jo 



dE st 



dE. { 



(33) 



The spectra of the pionic se condary e^, 7-rays, and neutri- 
nos are calculated using the iKamae et al.l d2006|) cross sec- 
tions f or proton energies below 500 TeV and the lKelner et aTl 
(2006) cross sections for proton energies above 500 TeV. This 
method has the advantage of being consistent between the dif- 
ferential cross sections and the pionic lifetime. We also con- 
sider other cross section parametrizations and pionic lifetimes 
in Appendix[A] we generally found that these produced simi- 
lar results. 

Calculating the injection rate of pair-produced e ± requires 
both the low-energy target photon and 7-ray spectra within the 
starburst. Assuming a planar geometry, with 7-rays traversing 
vertically out of the starburst disk 7 , the 7-ray number density 
is 



/V 7 (£ 7 ) 



Q-y(Ey)h 



[l-exp(-r 77 (£ 7 ))], 



(34) 



where Q 1 is the injection rate of 7-rays per unit volume in 
photons per unit energy per unit volume. For a photon num- 
ber density spectrum «(e), we calculate the pair production 
opacity t 77 using 



r 77 (£' 7 ) = / hn(e)a 11 (e,E 1 )de 



(35) 



where we use the approximation in lAharonianl d2004) (equa- 
tion 3.23) for the differentia l pair-production c ross s ection 
tr 77 (e,£ , 7 ) given in full in Gould & Schreder (1967]), and 
where we have also assumed that the radiation field is 
isotropic. When calculating the 7-ray spectrum observed at 
Earth, we replace h with £ ffi , the sightline to the center of the 
starburst disk (equal to R for perfectly edge-on disks). 

We use the GRASIL SEDs for the radiation fie lds of M82 
and Arp 220, scaled to the correct luminosities (Sil va et aTl 

7 Note that 7-rays do not always travel straight out of the disk plane, but 
also at horizontal angles through the disk. Thus we underestimate the mean 
r by a geometrical factor of order unity. 



[19981) . 8 For NGC 253, we use the M82 GRASIL SED scaled 
to the NGC 253 s tarburst's luminosity. We use the SED of 
iPorter etail d2008l) for the Galactic Center, adding an infrared 
greyb ody dust component fro m star-formation of temperature 
20 K dLaunhardt et al.ll2002l) . The GRASIL luminosities are 
converted to energy densities as U = L/(2irR 2 c). We then add 
the CMB to the SEDs of NGC 253, M8 2, and Arp 220' s nu- 
clei (the CMB is already present in the IPorter et al.ll2008l ra- 
diation field for the Galactic Center). We show the resultant 
7-ray optical depths for these SEDs in Figure Q] M82 and 
Arp 220 are both optically thick in the 10 - 100 TeV range. 
NGC 253's optical depth peaks at ~ 35 TeV with t « 0.5. 
However, the Galactic Center is transparent at all considered 
energies. Opacity from three radiation components are visible 
in FigureQ] near infrared radiation from old stars, far infrared 
radiation from dust grains, and the CMB. In NGC 253, M82, 
and Arp 220, the far infrared radiation completely dominates 
the 77 opacity. 

The source function 2pah(£e) of pair production <? ± 
is then calculated fr o m the 7-ray and IR spectra (e.g., 
lAharonian etail 119831: iBottcher & Schlickeiserl 1 1997b . We 
hav e implemented pair pro duction in our code. We use 
the lAharonian et al.l d 19831) approx imations for the source 
functio n, as given in equation 32 of Bottcher & Schlickeiserl 
d!997l) . which are known to be accurate at high energies. 

After calculating the 77 pair e ± spectrum, we then calcu- 
late the 7-ray luminosity from these pairs. We then in turn 
calculate the pair e ± from these 7-rays, using the same pro- 
cedure as we did for the 7-rays from protons and primary and 
secondary e ± . We found that these higher-order pair e ± were 
typically only a small fraction of the 77 pair e ± population, 
reaching a maximum of 50% in the Arp 220 nuclei models 
even with B = 250 fiG and under 8% for M82, NGC 253, and 
the Galactic Center. 

3.3. Constraints 

We run grids in 77, B, p, and the normalization of the diffu- 
sive escape time fdiff (see eq.l32i>. Hadronic models (with CR 
protons, secondary e , and the pair e associated with the 
7-rays from these) are run independently of leptonic models 
(with primary CR electrons and the pair e ± from the 7-rays 
they generate), giving us a hadronic and leptonic template for 
each parameter set. The two are then added together by scal- 
ing the hadronic template with rj and the leptonic template 
with £. For starbursts where there are error bars in the radio 
data (Galactic Center, NGC 253, and M82), we then select the 
free-free emission flux S t herm(GHz) at 1 GHz based on chi- 
square fitting of the radio data, scaling £ for each S± erm (GHz) 
so that the 1.4 GHz synchrotron radio emission of the star- 
bursts equals our mo del predictions. We use the interferomet- 
ric measurements in [Williams & Bower] d2010t) of M82 and 
the starburst core of NGC 253 an d the > GHz radio mea- 
surements in ICrocker et al.l d201 lbl) for the Galactic Center 
("HESS region") radio flux. We use the radio data compiled 

8 The GRASIL SEDs are ava ilable at 
http://adlibitum.oat.ts.astro.it/silva/grasil/modlib/fits/fits.html 

9 Note that both equations 26 and 32 in Bottcher & Schlickeiser 11997) 
are for electrons only (Aharonian & Khangulyan 20 1 1 , private communica- 
tion; Bottcher 2011, priva te communication). Our attem pts at calculating the 
pair spectrum using the Bottcher & Schlickeiser 1 1997) formula gives results 
that ar e a factor of 2 too small, both by comparing with the lAharonian et aT] 
1 1983) pair spectrum and checking energy c onserva t ion. A previous version 
of this paper assumed that Bottcher & Schlickeiser 1 1997) equation 32 was 
for both e ^ ; thus the pair injection rate was underestimated a factor of 2 in 
that version. 
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FIG. 1 . — On the left, we show the photon densities of the radiation fields of 
starburst nucleus of Arp 220 (dotted). The FIR emission from dust dominates t 
photons is comparable. On the right we show the pair-production optical depths 
M82 and Arp 220 are opaque at tens of TeV. 

in lTorresI ([2004) for the east and west nuclei of Arp 220; since 
these data do not have error bars, we simply normalize £ to 
match the observed 5 GHz flux. Models with negative £ (if 
the secondaries alone overproduced the radio flux) were not 
allowed. 

We then requi re the predicted 0.3 - 10 GeV jAbdo et alJ 
2010b) and TeV (lAcero et al.l 2009: Ac ciari et alJ2009h emis- 
sion of M82 and NGC 253 to match the observed values to 
within a factor 2. For the Galactic Center region, we require 
the predicted TeV em ission to match the observed emission 
dAharonian et alj|2006l) to within a factor 2, and the predicted 
GeV emission to be lower than the observed emission within 
\£\<l° and \b\ < 1° from EGRET dHunter et alJll997l) . For 
Arp 220, we require both the predicted GeV and TeV to be 
lowe r than the upper lim its from the F ermi-LAT one year cat- 
alog (lAbdo et alJl2010bl) and MAGIC (lAlbert et alJl2007h . re- 
spectively. 

We also run a purely calorimetric hadronic model for each 
p, B, and density, with no advective or diffusive losses. By 
integrating up the volumetric power generated in secondary 
e ± , 7-rays, and neutrinos in a model, and comparing the yield 
with the calorimetric model, we can quantify the efficiency of 
pionic losses as the calorimetric fraction F cs \\ 

"~ £°V<g. + E<g + E*g)dE- 

We also quantify the efficiency of pionic losses for generating 
VHE products by calculating F ca i(> 10 TeV), in which we 
change the lower bound of integration to 10 TeV. 

Usually, a range of r\ will be compatible with these con- 
straints. For the sake of brevity, for a given parameter set (E g , 
B, p, fdiff(3 GeV)), we consider only those models that either 
have (a) the minimum allowed r], (b) the maximum allowed 77, 
(c) the 77 which predicts the closest match to the ~ GeV 7-ray 
emission (in M82 and NGC 253) or the TeV 7-ray emission 
(in the Galactic Center), or (d) 77 fulfills (a), (b), or (c) for 
some other 7™™ and the model otherwise satisfies our con- 
straints. For Arp 220's nuclear starbursts, we simply assume 



he Galactic Center (short-dashed), NGC 253 (long-dashed), M82 (solid), and a 
le photon population, except in the Galactic Center where the number of CMB 
of these starbursts, along a vertical sightline out of the starburst disk. Note that 

that 77 = 0.1 since there are no 7-ray detections yet. 

4. RESULTS AND COMPARISON WITH OBSERVATIONS 

At GeV energies, synchrotron losses must compete with 
extremely strong bremsstrahlung, ionization, and IC losses, 
which cool e ± before they radiate much synchrotron. In or- 
der to account for the observed radio emission, there must be 
more CR e ± than one would naively expect for primary elec- 
trons with Milky Way-like acceleration efficiencies. In gen- 
eral, we find that two different kinds of models work, given 
the radio constraints and the 7-ray observations: 

1. ) In models with high B, the GeV CR e ± spectrum is 
dominated by pionic secondaries, which enhance the radio 
emission. In this limit, models are not affected much by vari- 
ations in sufficiently small £, since any small £ will result in 
sub-dominant primaries. The amount of synchrotron emission 
from secondaries (and tertiary pair e ± from pionic 7-rays) is 
set by the efficiency of CR proton acceleration (77 s» 0.1, in 
turn set by the 7-ray observations), the efficiency of pionic 
losses F ca \ and pair production (set by the gas density and ra- 
diation fields, not B), and the power of synchrotron with re- 
spect to other losses / sync h (which is highly dependent on B 
for GHz emission, but is ~ 1 for X-ray emitting energies in 
high B models). Thus the hadronic component asymptotes to 
a constant synchrotron X-ray luminosity as B increases. How- 
ever, the GHz radio synchrotron emission does increase with 
B, because of the dominant non-synchrotron losses; in models 
with too high B, the secondaries overproduce the radio emis- 
sion at fixed 77. This sets an upper bound on B. 

2. ) At low B, the synchrotron emission of secondary e ± 
is insufficient to explain the observed GHz radio. The only 
way for there to be enough CR electrons is to increase the pri- 
mary CR electron acceleration efficiency greatly. Thus, low B 
models favor primary electrons dominating the CR e ± spec- 
trum at GeV energies. 10 If the CR primary electron spectrum 

10 Lacki et al. 1 2010) postulated a "high-Eg conspiracy" that sets the ra- 
dio synchrotron luminosity of starbursts: the suppression of starbursts' radio 
emission by non-synchrotron losses is compensated by the appearance of sec- 
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ends at low energies, the 10 - 100 TeV CR e spectrum only 
consists of secondaries and tertiaries. If the primary CR elec- 
tron spectrum extends to higher energies, though, the primary 
electrons in these models may overwhelm the secondaries and 
tertiaries and greatly enhance the synchrotron X-ray emis- 
sion. Although we consider them, we feel that these models 
are contrived: in those starbursts with 7-ray constraints, these 
models require the CR proton acceleration efficiency to be the 
same or less than the in the Milky Way (77 < 0.1), so that the 
the total hadronic and leptonic 7-ray emission not exceed ob- 
servations. However, they also require that primary CR elec- 
trons be accelerated much more efficiently (£ w 0.05 - 0.2) 
in starbursts com pared to £ < 0.01 -0.02 in the Milky Way 
dLacki et al.ll2010t the numerical value depends on the shape 
of the CR electron spectrum at energies < GeV). They can 
be distinguished observationally from the high B secondary- 
dominant models by 7-ray emission below ^100 MeV: mod- 
els with low B have intense bremsstrahlung and IC emission 
which flattens out the pionic "bump" at lower energies. In- 
deed, the leptonic 7-ray emission sets a lower limit on B for 
the starbursts we model. 

This basic idea that different B require different elec- 
tron acceleration efficiency to match the GHz rad i o con - 
straints has been described befo re in iPersic et al.l d2008l) . 
Ide Cea del Pozo et al l d2009al) . and lRephaeli et all d2009t) 

In Table s [2j [3] EJ and[5] we list the synchrotron X-ray emis- 
sion in models for the Galactic Center, NGC 253, M82, and 
the nuclei of Arp 220, respectively. We find that the syn- 
chrotron X-ray emission is very weak in the Galactic Cen- 
ter and generally weak in NGC 253 and M82, especially in 
models where the primary e~ spectrum cuts off below a TeV 
(7m™ = 10 6 ). In Arp 220, however, the synchrotron contribu- 
tion to the diffuse X-ray emission is significant, even with a 
low energy primary cutoff. The IC emission is also a small 
minority of the X-ray emission. We discuss individual star- 
bursts in §[401-1441 

Some general trends are apparent in the tables. First, the 
synchrotron X-ray fraction is higher if the injection spectra 
are harder, simply because there are more e^ and protons at 
higher energy. Second, the amount of X-ray emission de- 
pends on the primary electron maximum energy. The cut- 
off dependence is especially strong in low B models where 
the primaries dominate the CR e spectrum. Even in high B 
models, where secondaries dominate the CR e ± population 
at ~ GeV energies, primaries can still dominate the ~ TeV 
e ± population, because diffusive escape becomes quicker at 
high energies and removes CR protons before they can in- 
teract with starburst gas. Third, a shorter diffusive escape 
time reduces the synchrotron X-ray contribution, particularly 
in high B models for M82 and NGC 253 where the secondary 
e ± at multi-TeV energies is dominant in some models and 
sub-dominant in others. 

To simplify our presentation, we also subjectively choose a 
"fiducial" model for each considered starburst. To qualify as 
a fiducial model, 8 must be near its approximate Milky Way 
value, within the range 50 - 100, whenever this is possible 
among the allowed models. This criterion selects the higher 

ondary e and the dependence of the critical synchrotron frequency on B. In 
the low B models, the non-synchrotron losses are still present and are even 
more influential because B is smaller, but they are instead compensated by 
a much larger primary CR electron acceleration efficiency than in the Milky 
Way. Thus, there still is a conspiracy, but a different one involving high £ 
instead of secondaries. 



B models. Because of the coarseness of our grids in B, how- 
ever, we were not able to get 5 to match in all of our fiducial 
models: 5 is 18 in the fiducial model for Arp 220's western 
nucleus, and 48 - 84 in the other starbursts. We also consider 
only p = 2.2 models (which also means 13 < S < 25, from 
eqn. I3TV Finally, we try to choose models that are close fits 
to any existing 7-ray observations. The parameters for these 
fiducial models are listed in Table Q] 

4.1. The Galactic Center 

Introduction - The Galactic Center region in many ways re- 
sembles a mini-starburst; we refer the reader to lCrocker et al.l 
d201 lbl) for a complete discussion of its properties. Dif- 
fuse hard (V « 2.3) TeV 7-rays have been observed with 
HESS in the region with \£\ < 0.8° and \b \ < 0.3° (R = 
112 pc ; h = 42 pc; lAharonian et all 120061: ICrocker et all 
1201 lbl) . which we take as our modelled region. The source 
of the CRs responsible for this emission has be en con- 
jectured to be Sgr A* ( e.g.. iBallantvne et al. 120071). super- 
novae in the region (e.g., Busching et al.ll2007t ICrocker et al.l 
1201 lbl) . a nd diffusive reacceleration of CRs in the inte rcloud 
medium (IWommer et all 120081: iMelia & Fatuzzoll201 ll) . The 
hard CR spectrum could either be ex plained by either a re- 
cent non-steady-state injection of CRs dAharonian et al.l2006l: 
iBusching et al.l 120071) . or s teady-state energy-in dependent 
transport such as advection (ICrocker et al.l |2Q 1 lbl). For the 
purposes of this paper, we assume like ICrocker et al.l (1201 lbl) 
that the diffuse TeV emission is from a steady-state popula- 
tion of CRs accelerated by star-formation processes, and then 
briefly discuss the consistency of this scenario with our re- 
sults. 

GeV emission from the nuclear starburst specifically, as op- 
posed to the inner Galaxy, is not yet well const rained, al- 
though it probably is present dCrocker et alj|20l"0ab . Radio 
emission is observed both from this region and a surrounding 
halo-like structure that extends out to \l\ < 3° and \b\ < 1° 
(R = 420 pc, h = 140 pc) which contains most of the observed 
radio emissi on. The luminosity of the dust in the region is 
4 x 10 8 L Q (lLaunhardt et al.ll2002l) : taking that as L TIR , we 
derive a star-formation rate of 0.07 M Q yr" 1 (eq. I2l1 i. which 
is compatible with th e star-formation rate over the past few 
million years found by Yusef-Zadeh et al. (2009) (see also the 
extensive discussion in ICrocker et al.ll201 1 alibi) . 11 This IR lu- 
minosity implies that the inner starburst region observed with 
HESS falls off the FIR-radio correlation observed in starbursts 
like M82 and N GC 253, but the sur rounding halo does lie on 
the correlation dCrocker et al.l20lTal) . This combined with the 
7-ray faintness of the starburst region given its supernova rate 
implies a strong wind, with a spe ed we take to be 600 km s" 1 
following Croc ker et al.l (1201 lbl) . The radiation energy den- 
sity from the dust luminosity is 42 eV cm" 3 . For the back- 
ground radiation field we add this star-formation TIR radia- 
tion fie ld to the inner Milky Way radiation field in Por ter et al.1 
(2008), which includes the contributions from stars outside 
the Galactic Center starburst proper. This radiation field has 
an energy density 15 eV cm -3 , for a total radiation energy 

11 lLaunhardt et al J <2002l) find a total stellar luminosity of 2.5 X 10 9 L Q , 
mostly from the nuclear star cluster, but this luminosity with eq. l21l implies a 
much higher star-formation rate of 0.3 Mq yr~' if it mostly came from young 
stars. Since the total radio emission is mostly from primaries in our models, 
this would require f to be ~ 6 times smaller than our derived values, or 0.001 . 
The TeV 7-ray emission, if it is hadronic, would imply ~ 6 times smaller rj 
than in our models, < 0.02. Crocker et al. (201 lb) however conclude that the 
supernova rate is consistent with the lower rate we use. 
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FIG. 2. — Our predicted total (synchrotron + f ree-free) radio spectru m (red, 
solid) compared to the observations compiled in Crocker et al. 1 201 lb]). Open 
circles are fitted while Xs are not fitted. The synchrotron spectrum itself is 
the red long-dashed line, with components from primary and secondary 
are shown as red dotted and red short-dashed lines. The black dash-dotted 
line is our fit to the thermal free-free emission with these parameters. 

density 57 eV cm" 3 (B rad = y/&irU md = 48 pG). 

The gas surface density that CRs traver se in the Galactic 
Center region is not very well constrained. iPierce-Price et al.l 
(2000) find a gas mass of 5 x 10 7 M(T) in the inner 200 pc , 
corresponding to E g = 0.09 g cm -2 . iFerriere et alj d2007l) 
find that S ? w 0.03 g cm" 2 in the inner few hundred pc, al- 
though steeply falling with distance from the Galactic Center. 
ICrocker et al.l d201 lbl) find that CRs have to sample gas of less 
than average density to fit the multiwavelength data. We run a 
grid in Y± g ranging from 0.003 g cm" 2 to 0. 1 g cm" 2 to account 
for these uncerta inties. We assume a magnetic field strength 
of 50 or 100 pG (ICrocker et alJl20T0ah . 

The Galactic Center region also has a ha rd X-ray and 
soft 7 -ray excess along the Galactic ridge (e.g.. lWorrall et al] 
1982). The central square degree has a 2 -10 keV sur- 
face brightness of 10~ 9 ergs s" 1 cm" 2 sr" 1 (Kova ma et al.l 
1996); we scale this to the modeled region, assuming 
constant surface brightness, to get a luminosity of 7.4 x 
10 36 ergs s" 1 . Although the X-ray ridge emission was once 
attributed to very hot gas, most of the emission has been re- 
solved into faint X-ray source s, mainly associated with stars 
dRevnivtsev et al]|2006l 12009). At higher energies > 50 keV 
(particularly so ft 7-rays), Invers e Compton may dominate 
the emission dPorter et al] 120081) . Synchrotron has previ- 
ously bee n considered as a source of the Galactic ridge X-ray 
emiss ion dProtheroe & Wolfendald[r980l: iPorter & Protherod 
1 1997b . but requires stronger magnetic fields than are present 
in much of the Milky W ay to avoid overproducing TeV In- 
verse Compton emission ( iAharonian & Atovanll2000|) . Given 
the strong magnetic fields and relatively high star-formation 
rate of the Galactic Center starburst, we reconsider whether 
synchrotron emission may contribute to the X-rays from this 
region. 

Modelled Radio and 7-rays - The Galactic Center region 
has few radio data points, so it is not surprising that we are 
able to fit the > GHz radio data. In the allow ed models, the 
radio emission is dominated by primaries (c.f. ICrocker et al] 



201 lb), which have £ = 0.005-0.01 for B = 100 //G and f = 
0.03-0.06 for B = 50 fiG. The 74 and 330 MHz data points 
fall well below our fits, which is probably due to the strong 
free-free absorptio n known to be present in the region (e.g., 
iBrogan et al.ll2003l) . 

Models with S g = 0.003 g cm" 2 are strongly preferred by 
the 7-ray constraints (see Table [2] for the allowed mod els). 
This supports recent modeling bv lCrocker et al] d201 lbl) im- 
plying that CRs do not sample the average density in the 
Galactic Center region. Low gas surface densities means that 
F cdi is also small, ~ 0.3-0.8% for E ? = 0.003 g cm" 2 (contrast 
this with f ca i sa 8-21% for E ? = 0.1 g cm" 2 ). Thus, in many 
of our models, the Galactic Center is actually less of a proton 
calorimeter than the Milky Way as a whole. Above 10 TeV, 
diffusive escape can reduce this even further, to ~ 0.4% when 
f di ff(3 GeV) = 10 Myr and - 0.06% when f diff (3 GeV) = 1 Myr 
for E g = 0.003 g cm" 2 . Secondaries are therefore a minority 
compared to the primaries (dashed lines in Figure[5J left). Ter- 
tiary pair e ± are negligible at all energies, because the Galac- 
tic Center is transparent to 10 - 100 TeV 7-rays (Figure Q] 
short-dashed line). 

Synchrotron X-rays - Except in a few of the most opti- 
mistic models, only a small fraction of the hard X-ray emis- 
sion from the Galactic Center is synchrotron (see Table [2]). 
This is perhaps not surprising, given that most of the X-ray 
emission has already been resolved. With CR protons and 
VHE 7-rays escaping so readily, there are few secondary or 
tertiary e at 10 - 100 TeV energies to emit synchrotron X- 
rays (Figure[3]l. Therefore the X-ray emission strength is de- 
termined by primary electrons. Furthermore, the relatively 
low B of the Galactic Center compared to the other star- 
bursts means that a given synchrotron frequency corresponds 
to higher e ± energies, where the cutoff in CR protons causes 
a steep falloff in the secondary e spectrum. We find syn- 
chrotron fractions of 0.01% to 6%, with a fiducial value of 
5.9% (Lf^ = 4.3 x 10 35 ergs s" 1 w 3 x 10" 7 Dm)- This is 
consistent with the resolution of 88 ± 12% of the X-ray ridge 
emission at energies of ~ 7 keV dRevnivtsev et al.l[2009h . For 
p = 2.2, we find 7m™ = 10 9 models have synchrotron frac- 
tions of ~ 6%; almost all allowed models with 7 m ™ = 10 6 
have synchrotron fractions of 0.1% or less. 

The synchrotron emission has a very soft spectrum (r2-io = 
2.35-2.50) in models with 7m™ = 10 6 . This is because in the 
weaker 50- 100 pG magnetic fields of the Galactic Center, the 
2 - 10 keV X-ray emission traces e of very high energy (grey 
shading in left panel of Figure[3]). The CR secondary e ± spec- 
trum starts to cutoff at these energies, because the CR proton 
spectrum ends at a PeV in our model. In reality, the "knee" 
of the CR proton spectrum may be at a different energy in the 
Galactic Center region, so the softening of the synchrotron 
spectrum may not be real (see Appendix [B] where we find 
harder hadronic X-ray spectra if 7 max = 10 7 ). In models with 
7m™ = 10 9 , the synchrotron spectrum has IVio = 2.12-2.27, 
since the primary electron spectrum remains hard at these en- 
ergies (left panel of Figure [3). 

Inverse Compton X-rays - The IC contribution to the ob- 
served 2-10 keV emission from the TeV-emitting region is 
also small. The main factor affecting the IC contribution is 
B: in models with B = 50 pG, IC makes up 0.5% of the 2 - 
10 keV emission; in models with B = 100 fiG, IC makes up 
0.09 - 0.16% of the 2-10 keV emission. The IC emission 
has a spectral index of 1 .45 - 1 .67, due to the strong advective 
losses flattening the electron spectrum. 
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FIG. 3. — The CR electron (left) and photon (right) spectra of the Galactic Center region. On the left, primary electrons are solid (7m™ = 10 6 ) and long-dashed 
(Tmax 1 = 10 9 ), pionic secondaries are short-dashed (blue is e + , black is e~, and grey is total). Pair and knock-off electrons are insignificant. Blue and black shading 
denote the approximate e*- energies that radiate in the 2-10 keV band through IC and synchrotron, respectively. The vertical line indicates the approximate 
energy that radiates synchrotron at 1 .4 GHz. On the right, synchrotron is red, bremsstrahlung is green, IC is blue, pionic is violet, and total is black. Solid lines 

are the Earth- observed absorbed luminosities in the 7 m '™ = 10 6 models. Dashed lines show 77 absorbed luminosities in Tniax = 10 9 models. The TeV data points 
are from Aharonian et al. 12006). 



The ratio of synchrotron and IC emission in the 2-10 keV 
band has a large range (0.029 - 55; fiducial: 48) in the al- 
lowed models. However, in all of the 7^™ = 10 9 models, syn- 
chrotron is brighter than IC, with typical ratios of ~ 3-5 for 
p = 2.4 models and ~ 50 for p = 2.2 models. The large values 
of these ratios even when p = 2.4 arises because Ub > f/ ra d 
and winds remove electrons below 10 GeV before they can 
radiate much IC (eqns. [9] and 1121 . By contrast, in most of 
the 7m 1 ™ = 10 6 models, the synchrotron to IC ratio is almost 
always less than unity. 

Conclusion - The Galactic Center in some respects is a poor 
target for synchrotron X-ray emission. This is because it has a 
large X-ray background probably from stellar X-ray sources, 
contains a fast wind that removes protons before they can in- 
teract with the gas and produce secondary <? ± and 7-rays that 
could cascade, and a weak radiation field compared to other 
starbursts that means few pairs are produced from the 7-rays. 
The synchrotron emission is a few percent of the observed 
emission only if the primary electron spectrum extends to high 
enough energies. IC emission in the 2-10 keV X-ray band 
is also very weak (< 1%), because the Galactic Center wind 
removes GeV electrons before they can IC upscatter ambi- 
ent radiation. However, the proximity of the Galactic Center 
means that many point sources can be resolved out, leavin g 
the truly diffuse X-ray emission (e.g.. lRevnivtsev et al.12 009). 

A potential concern about our results is whether time- 
independence is a good approximation. Our assumed TIR 
luminosity of 4 x 10 s L Q translates to a supernova rate of 
(2200 yr)" 1 , while our fiducial diffusive escape time is fdiff = 
10 6 yr (E/3 GeV)" 1 / 2 = 1700 yr (E/PeV)" 1 / 2 . Thus the 
steady-state approximation roughly holds for protons in our 
models until a PeV. The advective lifetime in our models, 
63000 yr (c.f. eqn.[T2"ii is long enough for ~ 30 supernovae to 



go off in the HESS region, which would smooth out stochastic 
effects at energies where advection dominates. Advective es- 
cape al so naturally produces the hard spectrum (ICrocker et al.l 
1201 lbl) . We cannot rule out that diffusive escape is much 
faster than we suppose, though. For primary electrons, the to- 
tal > 10 TeV population of the entire region is not affected by 
diffusive escape since the synchrotron lifetime alone at these 
energies is f synch = 120 yr^/lOTeVrHB/lOO^G)- 2 (eqn.[3), 
faster than any possible escape process. However, the level of 
primary CR electrons may vary overkyr timescales if the time 
in which CR accelerators accelerate 10 - 100 TeV electrons is 
shorter than ~ 2200 years. Thus, there could be a stochastic 
element in the synchrotron X-ray emission from the Galactic 
Center region. 

Another issue with these is that the Galactic Center is em- 
bedded in a 420 pc wide radio halo, in turn embedded in the 
Galaxy as a whole. It is therefore possible that our one-zone 
approximations cannot capture the complex geometry of the 
region. We generally expect the background sea of CRs from 
the Galaxy to be only a minor part of the CR p opulation in the 
R < 1 12 pc region. As I Aharonian et ail (j2006) noted, the TeV 
detection indicates a CR energy density several times greater 
than that of the Galaxy at large (~ 1 eV cm 3 ). Using t h eir own 
one-zone models of the same region, ICrocker et alJ JMTTbh 
also find that CR energy densities of 2.5-50 eV cm" 3 are con- 
sistent with the multiwavelength data, with the best-fit model 
having Ucr = 20 eV cm -3 . 

In the models presented here that are consistent with con- 
straints, we find CR energy densities of 1.8-38 eV cm" 3 . 
In all of these models, then, the Center-accelerated CR pop- 
ulation is at least comparable, and often overwhelms any 
background Galactic population. The models at the low end 
of the range have small 77 (~ 0.01 compared to the canon- 
ical 0.1) and TeV luminosities that are relatively small. If 
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we restrict our attention to only those models that also have 
0.05 <r/< 0.2 (within a factor of 2 of the typically assumed 
0.1), we find CR energy densities of 4.9-37 eV cm -3 . Our 
fiducial mode l has a C R energy density of 20 eV cm" 3 , just as 
ICrocker et all (1201 lbl) find. 

Since we are interested in the high energy electron pop- 
ulation, we also compare the energy density in CRs with 
> 10 GeV. In the local Mil ky Way, this en ergy density is 
just 6 x 10" 4 eV cm" 3 (see lAdriani et al.ll201 ll and references 
therein). In the models consistent with constraints, however, 
the > 10 GeV e ± energy density is 0.016-0.10 eV cm" 3 , 
with a value of 0.035 eV cm" 3 in our fiducial model. These 
high energy densities occur even in the face of the strong syn- 
chrotron and IC losses that CRs experience, requiring a large 
luminosity density of CR e to sustain. We therefore con- 
clude that the Galactic Center CR population is distinct from 
that of the large-scale Galaxy. 

The presence of the 420 pc radius radio halo which actu- 
ally has a larger GHz synchrotron luminosity than the central 
molecular zo ne may also pose a cha llenge for our use of one- 
zone models. ICrocker et alj d201 lbl) interpret this halo as be- 
ing the emission from CR e advecte d out of the small e r TeV- 
emitting region within. In addition, Launhardt et al. (2002) 
find this halo region emits an additional 2 x 10 8 L Q of bolo- 
metric luminosity, which translates to a star-formation rate 
half that of the inner 112 pc. Multi-TeV e generated in 
the inner region do not live long enough to contribute to any 
synchrotron X-ray luminosity from the radio halo. However, 
GeV e ± that can emit IC X-rays are likely present. To esti- 
mate how big of a foreground IC X-rays from the radio halo 
is, we find the ratio of 1 .4 GHz radio surface brightnesses of 
the halo region (\£\ < 3°, \b\ < 1°) a nd the central molecu- 
lar zone (\£\ < 0.8°, \b\ < 0.3°) from ICrocker etail (1201 lbh . 
The central molecular zone has a 1 .4 GHz surface brightness 
that is ~ 3 times that of the outer radio zone. Thus, while the 
central molecular zone CR population dominates the observed 
radio emission in that area of the sky, the radio halo also con- 
tributes significantly. A multidimensional model of the region 
is necessary to fully capture the CR transport in the region, 
although for order-of-magnitude purposes our one-zone ap- 
proach is sufficient. 

4.2. NGC 253 Starburst Core 

Introduction - NGC 253 contains one of the nearest star- 
bursts (we adopt D = 3.5 Mpc; 1 ."= 1 7 pc) and one of the two 
detected in both GeV (by Fermi-LAT; lAbdo et al.ll2010at) and 
TeV 7-rays (by HESS: lAcero et all2 009). Previous modelling 
of the nonthermal em i ssion from NGC 253 has been done by 
Pa glio ne" et al.1 ([1996]) . iDomin go- San tamarfa & Torresld2005l) . 
and iRephaeli et al.l d2009l) . and the latter predictions are in 
good agreement with the Fermi and HESS results. Unlike 
the Galactic Center, the starburst core of NGC 253 lies on 
the FIR-radio correlation and is 7-ray bright compared to the 
Milky Way. NGC 253 has a large reservoir of gas in its cen- 
ter; the estima ted mass is (2-5) x 10 7 Mq from molecu- 
lar lines (e. g., Mauersber ger et ail 1996tlHarrison et al.ll 19991: 
iBradford et al.1 120031; ISakamoto et al.l 1201 ll) . corresponding 
to central gas surface densities of 0.07-0.15 g cm" 2 . The 
starburst appears as a disk with a di ameter of 20"- 30 
" (340 - 510 pc) in radio continuum dTurner & Hoi 119831; 
lUlvestad & Antonuccil fl997l) an d 15"- 30" (250 - 510 pc) 
in inf r ared and molecular lin es dPeng et al.1 119961: lUlvestad! 
l2000t ISakamoto etail 1201 lb ; we adopt a radius of 150 
pc. The vertical extent of the radio disk is 8" (136 pc) 



dUlvestad & Antonuccill 19971) . implying a maximum vertical 
scale height of 68 pc. On the other hand, a cylindrical geome- 
try for NGC 253's radio disk, given an inclination of 78° and 
a radio diameter of 340 pc, implies that ^ 7 1 pc of the vertical 
extent is the projection of the diameter; the remainder implies 
a scale height of only 33 pc. Due to the uncertainties in the ge- 
ometry, we adopt a scale height of 50 pc (giving a gas density 
of 190 cm" 3 ) Finally, the total IR luminosity of NGC 253 is 
4 x 10 10 Lq (San ders et al"1l2003l). of w hich half comes from 
the nuclear starburst (Melo et al. 2002). We therefore adopt a 
starburst bolometric luminosity of 2 x 10 10 Lq, which corre- 
sponds to a low supernova rate of 0.025 yr" 1 , and a radiation 
energy density 1200 eV cm" 3 , equal to that of a magnetic field 
B md = 220 nG. 

NGC 253 ha s also been studied with Chandra and XMM- 
Newton (e.g.. JStrickland et al l 120001: IWeaver et al.1 120021: 
IStrickland et al.l2002l:lBauer et al.l2008l) . revealing the effects 
of a starburst wind in soft X-rays and hard emission from 
the galactic disk. XMM-Newton measured the diffuse (point 
source subtracted) 2 - 10 keV luminosity from the optical disk 
as 8.5 x 10 38 ergs s" 1 dBauer et al.ll2008l) . We adopt this as the 
X-ray luminosity to compare against, although it includes the 
regions outside the starburst ; these outer regio ns host 30% of 
the total star-formation rate (iMelo et al.ll2002l) . 

Modelled Radio and 'y-rays - The combination of radio, 
GeV, and TeV constraints still leaves a large and complex 
parameter space of models that work. Selected models are 
shown in Table [3] 

We show radio emission in a typical model in Figure@] The 
flat radio spectrum favors high GHz thermal fractions in our 
models, with a span of ~ 4-28% (fiducial: 21%) selected 
by our criteria. As seen in Figure [4] there is evidence for 
unaccounted spectral curvature in the residuals, and the high 
thermal fraction overproduces the observed emission at 20 - 
100 GHz. This implies that the synchrotron spectrum is actu- 
ally flatter than predicted by our models at 1 GHz, and then 
grows steeper at higher frequency, with a lower thermal frac- 
tion. The flatter spectrum can arise from loss mechanisms 
with a flatter energy dependence than synchrotron or IC, such 
as bremsstrahlung, ionization, or advection losses. Enhancing 
advection losses decreases the yield of secondary electrons, 
but increasing the density (and bremsstrahlung and ioniza- 
tion losses) increases the yield. Furthermore, more primaries 
would be needed to account for the smaller fraction of elec- 
tron power going into GHz synchrotron emission if there are 
additional losses, so the results of such changes are complex. 
Secondaries dominate the GHz synchrotron emission in mod- 
els with high B, but in low to intermediate B, primaries domi- 
nate. 

The high energy emission from example models is shown 
in Figure [5] The predicted GeV 7-ray emission changes from 
being mainly leptonic at B = 50 fiG to nearly all hadronic 
at B = 150 /iG. This implies that NGC 253's starburst can- 
not have a magnetic field much lower than 50 /iG, or else 
the leptonic 7-ray emission would exceed the observed val- 
ues. Similarly, the magnetic field cannot be much higher than 
150 /iG, given our assumed density and radiation field, or else 
the radio emission from the secondaries would be overpre- 
dicted, and in fact no models with B > 200 /iG fit our crite- 
ria. Similar values have been derived from previous modeling 
of NGC 253's starbu rst dDomingo-Santamarfa & Torresl2005l 
IRephaeli et al"1l2009l) . 

In B > 100 /iG models, high 77 models are favored, and our 
fiducial model has r\ = 0.40, three to four times higher than 
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FIG. 4. — Our predicted total radio spectrum for the sta rburst core of NGC 
253 (r ed, solid) compared to observations compiled in Williams & Bower 
(2010) (open circles). The line styles are the same as in Figure [2] The ra- 
dio emission at high frequencies is overproduced by the large thermal com- 
ponent, suggesting that the synchrotron spectrum is intrinsically flatter than 
shown here and that escape, bremsstrahlung, or ionization are stronger and £ 
is higher in reality. 



the other starbursts. This is because the proton calorimetry 
fraction is less than one half of the value for M82 in our fidu- 
cial model, even though the 7-ray to bolometric luminosity 
ratio is the same. The actual CR acceleration efficiency could 
be ~ 0.1 and still consistent with the 7-ray data if: (1) the 
supernova rate in NGC 253's core is higher than we assume 
(and 5 is the same as in our fiducial model), since only the 
injection rate of CRs matter, (2) the gas density is higher than 
we assume, increasing F ca i, or (3) adyective losses are weaker, 
again increasing F ca i (c.f. lLacki et al.ll20TTb . 

The HESS detection also constrains the proton and electron 
populations at TeV energies. In low B leptonic models, the 



primary electron spectrum cannot extend to 7 ma x =10 with 
p = 2.0, or else the IC emission is overproduced. Even with 
p = 2.2, the leptonic models are severely constrained and just 
barely fit the TeV constraint (see Table [3). The high B sce- 
narios allow a few p = 2.0 models, but only when diffusive 
escape time is quick and the TeV/GeV ratio is still ~ 3 times 
higher than observed. 

In our models, NGC 253 is not a true proton calorimeter be- 
cause of its strong wind, but it is much more calorimetric than 
the Milky Way. The calorimetry fraction obtained by inte- 
gratin g the pionic pro ducts over all energies is F ca i w 10— 30% 
(c.f jLacki et al.l201 ll) . However, energy-dependent diffusive 
escape plays an increasing role at higher energies. Above 10 
TeV, F ca i drops to ~ 12% when f diff (3 GeV) = 10 Myr and 
~ 2% when f diff (3 GeV) = 1 Myr. 

Synchrotron X-rays - The wide span of the parameter space 
allowed by the data (Tabled means the 10 - 100 TeV electron 
population is not very well constrained, even with the HESS 
detection. Thus, the synchrotron X-ray fraction of NGC 
253's diffuse X-ray emission ranges from 0.19% to 120%, 
with a fiducial value of 8% (L^ = 7.2 x 10 37 ergs s" 1 
9 x 10~ 7 Ltir). For p = 2.2 models, the synchrotron fraction 
is typically tens of percent in low B models with high 7m™, a 



few percent in high B models, and a few percent to less than 

one percent in low B models with low 7m™. Models with 

large synchrotron fractions (> 1/3) overproduce the observed 

TeV emission by a factor ~ 2. 

The spectral index of synchrotron emission is 2.28 < 

T 2 -io < 2.42 for 7S"™ = 10 6 and 2.00 < r 2 _i < 2.29 for 
prim _ , 9 
/max — lu . 

Inverse Compton X-rays - The calculated IC X-ray emis- 
sion in the 2-10 keV band varies because the magnetic 
field in NGC 253's starburst is unknown. The models with 
B = 50 /iG, 100 pG, and 150 ^G have IC fractions of 7- 19%, 
2-4%, and 1.3-1.9%, respectively. The synchrotron-to-IC 
ratio varies widely from 0.013 to 49 (fiducial: 2.8). Mod- 
els with high 7™™ have higher synchrotron/IC ratios (see Ta- 
ble 0. The IC emission is consistently much harder than the 
synchrotron emission with 1 .26 < T2-io < 1 -40. 

Conclusion - The GeV and TeV detections of NGC 253's 
starburst inform our knowledge of the synchrotron X-ray 
emitting e ± population, but leave room for large variations. 
While the synchrotron fraction of NGC 253 can be as low 
as 0.2% or as high as 120%, in the fiducial model it is 8%. 
A key uncertainty is the primary electron energy cutoff, be- 
cause the <? ± population is dominated by primary electrons at 
these energies. The IC emission is also a minority of the X- 
ray emission, with the largest uncertainty being the unknown 
magnetic field in NGC 253. 

The higher resolution measurements of NGC 253 with 
Chandra can help reduce the competing emission from other 
sources. First, Chandra can resolve out more point sources 
with its higher angular resolution. Second, the emission from 
the outlying disk can be subtracted, leaving only the emis- 
sion from the nuclear starburst proper. Chandra has revealed 
that there is extended hard X-ray emission in the nuclear star- 
burst of NGC 253, but litt le is known about it (IStrickland et alJ 
l2lMIWeaverltini2l)02h . 

4.3. M82 

Introduction - M82 (D = 3.6 Mpc; 1"= 17 pc) is also de- 
tected in GeV and TeV 7-ray s , with Fermi-LAT and VER - 
ITAS respectively (lAbdo 'et all I20T0H lAcciari et alJ l2009h . 
Ide Cea del Pozo et all d2009al) and iPersic et alJ d2008h pre- 
viously modelled the nonthermal emission from M82, and 
there is again good agreem ent with the observed fluxes 
Ide Cea del Pozo et al.1 ©09b). M82 lies on the FIR-radio 
correlation and is 7-ray bright with respect to the Milky Way. 
Most of the star-formati on is in the starburst in the center of 
M82. IWeifi et al.1 (1200 ll) find a gas mass of 2.0 x 10 8 M Q in 
the inner 280 pc (for D = 3.6 Mpc), corresponding to a gas sur- 
face density of 0.17 g cm" 2 , which we use here. The infrared 
emission of M82 is concentrated in a ring with radius 225 
pc, and the gas p eaks at a radius of 250 pc (see the compila - 
tion in Table 4 oflGoetz et al.| [l990: see also lKennicuttlll998h . 
I Williams & Bower! d2010l)~ find a radio extent for M82 of 35." 
by 10", corresponding to a radius of 300 pc and a scale height 
of 90 pc. We therefore adopt 250 pc as the radius and 100 
pc as the scale height of the modeled starburst re gion. Other 
studies find radio scale h eights of 50 - 200 pc (Kle in et alJ 
119881: ISeaauist et al.ll 19851) . Values for the wind speed vary; 
however, iGrevd d2004l) find that the win d accelerates over a 
200 p c scale height to 400 km s" 1 and Westmoauett e et all 
(2009) find Ha FWHM of 100-300 km s" 1 in the inner 
few hundred parsecs. The wind speed may asymptote to 
> 1000 km s" 1 , but in the standard theory this happens out- 
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FIG. 5.— The CR electron (left) and photon (right) spectra of NGC 253's core for models with p = 2.2, f dif f(3 GeV) = 10 s yr, and 77 = 0.40. On the left, 
primary electrons are solid (7m™ = 10 6 ) and long-dashed (7m™ = 10 9 ), pionic secondaries are short-dashed (blue is e + , black is e~, and grey is total), knock-off 
electrons are dotted, and pair-production tertiaries (using the 7-ray spectrum for 7m™ = 10 9 ) are dash-dotted. Note that pair production are comparable to 
pionic secondary at 10 TeV. Blue and black shading denote the approximate energies that radiate in the 2-10 keV band through IC and synchrotron, 
respectively. The vertical line indicates the approximate energy that radiates synchrotron at 1 .4 GHz. On the right, synchrotron is red, bremsstrahlung is green, 
IC is blue, pionic is violet, and black is total. The dotted lines are without 77 absorption (7 m '™ = 10 6 ), while solid are 77 absorbed according to the sightline 

observed from Earth (7m™ = 10 s ). Dashed lines show 77 absorbed luminosities in 7m™ = 10 9 models. Note that the B = 50 fiG model with 7m™ = 10* is not 
allowed by our constraints. Observed luminosities from from Fermi and HESS (triangles), as well as the diffuse emission from XMM-Newton (square) and total 
emission from BeppoSAX (star; Cappi et al. 1999), are plotted. X-ray luminosities are scaled to one In bin in energy assuming T = 2.0. The effects of neutral 
hydrogen absorption are not shown on right. 



side most of the starburst; in side the starburst we expe ct wind 
speeds near the sound speed ([Chevalier & Clegg 1985), which 
is ~ 200-300 km s" 1 dGrevd200 4). We therefore again adopt 
an advection speed of 300 km s" 1 . The total IR luminosity of 
M82 is 5.9 x 10 10 L from lSanders et alJ (120031) . correspond- 
ing to a supernova rate of 0.06 yr" 1 . This also implies a radia- 
tion energy density 1300 eV cm" 3 , equal to that of a magnetic 
field B mA = 230 ^G. 

The diffuse X-ray emission of M82 is well-studied. The 
diffuse hard emission has a luminosity of 4.4 x 10 39 ergs s" 1 
and is spatially extended over a few hundred pc, with a scal e 
height estimated at h = 175 pc (Stricklan d & Heckmanl2007h . 
The spectral slope of the diffuse hard X-ray emission is T« 
2-3 and it is difficult to explain as therm al bremsstrahlung or 
IC emission (IStrickland & Heckm an 2007). 

Modelled Radio and j-rays - The GHz synchrotron ra- 
dio spectrum combined with the GeV and TeV 7-ray detec- 
tions leaves us with the parameter space in Table [4] The 
predicted synchrotron radio emission tends to be too large 
below a GHz (as s hown in Fig ure H); this may be due to 
free-free absorption (Kle in et al.lll988l) . Conversely, the syn- 
chrotron radio spectrum falls off steeply at high frequencies, 
with a » 0.8-0.9 instead of the observed 0.7 ( Kl ein et all 
ll988t fWilliams & Bower 2010). In order to accommodate the 
observed spectrum, our fitting requires higher 1 GHz ther- 
mal fracti ons (8 - 21%; with a fiduc ial value of 15%) than 
derived in IWilliams & Bowerl (1201 Ol) (6 ±2%). These high 
thermal fractions overproduce the observed total radio emis- 
sion at 20 - 100 GHz. The underlying synchrotron spec- 
trum is therefore probably flatter, which could occur with 
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FIG. 6. — Our predicted sync hrotron radio spectrum (re d, solid) com- 
pared to observations compiled in Williams & Bower (2010) (open circles) 
and the unfit data in Klein et al. 1 1988) (Xs). The residuals indicate the syn- 
chrotron spectrum is too steep at high frequencies, suggesting that escape, 
bremsstrahlung, or ionization are stronger and £ is higher in reality. The line 
styles are the same as in Figure[2] 

higher bremsstrahlung, ionization, or escape losses. When 
B > 150 /iG, the GeV e ± spectrum is dominated by pionic 
secondaries (short-dashed lines in left panel of Figure |7), and 
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the > 100 MeV emission is almost entirely pionic 7-rays (pur- 
ple line in right panel of Figure[7]l. These models have low £, 
and thus relatively little power going into primary electrons. 
When B < 100 fxG, however, the GeV e ± spectrum is mostly 
primary electrons (long-dashed line in Figure [7]); these mod- 
els usually have much higher £ (6 < 10 for p = 2.2) than the 
Milky Way. 

As with NGC 253, the GeV 7-ray emission is primarily lep- 
tonic when B = 50 /iG and hadronic when B = 150-200 /iG. 
The predicted GeV 7-ray flux is already at least 40% higher 
than observed in the B = 50 /iG models because of the strong 
bremsstrahlung and IC emission; therefore, M82's magnetic 
field is probably greater than 50 /iG. Indeed, most of these 
models only work by either cutting off the primary <? ± spec- 
trum so that the IC emission doesn't greatly overproduce the 
observed TeV emission, or by having very low r\ so there's 
virtually no pionic emission and all of the TeV emission is 
leptonic. Thus, the TeV detection of M82 place constraints on 
the high energy electron spectrum when B is low. Conversely, 
B is unlikely to be much higher than 200 /iG (if our assump- 
tions about the gas density and radiation field are correct), or 
else the secondary e that accompany the pionic 7-rays pro- 
duce too much radio emission. These values for B are in ac- 
cord with previous modeling for M 82's starburst dPersic et alJ 
I2008t Ide Cea del Pozo et al.ll2009ah . The efficiency of proton 
acceleration can be small in the low B leptonic models, but 
when B > 100 /iG and p = 2.2, 7/ is of order 0.1, the usual 
value assumed in models of star-forming galaxies. 

In our models, M82 is slightly more proton calorimetric 
than NGC 253, with F cal sa 10-50 % integrated over all ener- 
gies. At the energies above 10 TeV relevant for synchrotron 
X-ray production, F ca \ is only ^11% when fdiff (3 GeV) = 
10 Myr and ~ 1% when f diff (3 GeV) = 1 Myr. 

Synchrotron X-rays - In Table @] we list the predicted X- 
ray synchrotron emission for a variety of model parameters. 
Overall, the results are similar to those for NGC 253: the dif- 
fuse synchrotron fraction varies from 0.4% to 54% but has a 
fiducial value of ~ 2% (Li^f 1 = 9.7 x 10 37 ergs s 1 ; assuming 

T w 2.0, Lj-To 1 ~ 5 x 10 -7 Ltir)- For p = 2.2 models, the syn- 
chrotron fraction is typically ~ 15% for low B models with 
a high energy primary e^ cutoff, 0.8-5% for high B mod- 
els, and a fraction of 1 % for low B models with a low energy 
primary e^ cutoff. Like NGC 253, most of the models at the 
high end of the range of synchrotron fraction have TeV lumi- 
nosities that are > 1 .5 times higher than observed. 

Interestingly enough, the spectral slope of synchrotron 
emission in the 2-8 keV energy range (V i=a 1.98-2.34) also 
matches the observed diffuse hard X-ray emission (V sa 2-3). 

Inverse Compton X-rays - The IC fraction of the diffuse X- 
ray emission is largest (4 - 10%) in the B = 50 /iG models and 
becomes smaller as B increases: 1.4 - 2.4% when B = 100 /iG, 
0.8-1.1% when B = 150 pG 0.6% when B = 200 /iG This 
confirms the expectation in Stricklan d & Heckmanl (120071) 
that IC does not make up the observed diffuse hard X-ray 
emission in M82. Because of the poorly constrained syn- 
chrotron X-ray fraction, the synchrotron-to-IC ratio is any- 
where between 0.05 and 36, with a fiducial value of 2.3. We 
find that the IC emission has 1 .23 < T2-8 < 1 -36, significantly 
harder than the observed spectral slope of the diffuse hard 
emission in M82. 

Conclusion - The synchrotron fraction of M82's diffuse 
hard X-ray emission is largely unknown, even with the GeV 
and TeV detections, but is of order 2% in our fiducial model. 



The unknown primary electron energy cutoff is one of the 
largest contributors to the uncertainty. IC is a minority 
(0.5 - 10%) of the diffuse hard X-ray emission, with its con- 
tribution depending mainly on magnetic field strength. The 
synchrotron-to-IC ratio in our models is ~ 2 in our fiducial 
model. 

Unlike NGC 253, foreground hard X-ray emission from the 
host galaxy has already been removed, and the extant stud- 
ies with Chandra already take into account the bright point 
sources. However, the diffuse hard X-ray emission has a dif- 
ferent morphology than the radio and IR disks (and therefore, 
the expected sources of CRs) with h sa 175 pc. Restricting 
attention to the radio disk would help studies of the nonther- 
mal contribution. Furthermore, the IC contribution might be 
identified spectrally since it has a hard spectrum. Finally, the 
models with the largest synchrotron X-ray contributions are 
dominated by primaries; so large synchrotron fractions could 
be tested by looking for concentrations of X-ray emission near 
CR accelerators. 

4.4. Arp 220 

Introduction - Arp 220 consists of a starburst disk con- 
taining most of the gas and possibly most of the star- 
formation, and two inten se but smaller starburst nuclei 
dDownes & Solomon] 119981) . While there is radio data for 
the individual starbursts, and the galaxy as a whole, there are 
no flux values in the literature for the starburst disk specifi- 
cally. We therefore consider the nuclear starbursts separately 
and ig nore the surroundin g disk. We adopt D = 79.9 Mpc 
from [Sanders et al. (2003]). The bolometric lu minosity of 
Arp 220 is 1.6 x 10 12 L (ISanders etalJl2003h : if it is en- 
tirely due to star-formation, this would imply a supernova 
rate of 2 yr" 1 . Estimates of the bolometric luminosities of 
the nuclear star bursts vary, but we a dopt 3 x 10 11 L Q for 
each , based off [Downes & Solomon! (Q998) and consistent 
with lSakamoto et all d2008l) . These imply supernova rates of 
0.34 yr" 1 in each nucleus. We note that some authors argue 
that the bolometric luminosity of Arp 220's western starburst 
is dominated by an AGN, based on its compactness, which 
would reduce the IR lum inosity attributable to star-formation 
(D ownes & E ckart 2007). On the other hand, the actual su- 
pernova rates may in fact be up to near 3 yr" 1 in the west- 
ern starburst and 1 yr" 1 in the eastern s tarburst, based on the 
appeara nce of radio-brigh t supernovae (Lon sdale et al.ll2006t : 
earlier. | Rovilqs et alJ 120051 found 0.7 yr" 1 ). According to 
lLonsdale et alJ d2006l) . the radio sources in Arp 220's nuclear 
starbursts are concentrated into regions ~ 100 pc wide, and 
the radio images of the n uclear starbursts indica te diameters 
of 80- 120 pc. Following lSakamoto et all fl999). we assume 
that the nuclear starbursts are small disks, each with radius 
and scale height 50 pc. The radiation energy density in each 
nuclear starburst is then 1.6 x 10 5 eV cm" 3 (B ra d = 2.5 mG). 
The gas mass in each disk is ~ 10 9 within a 100 pc radius 
of each dDownes & Solomonll998l:ISakamoto et aljl999l) . for 
gas surface densities of 7 g cm" 2 . Given the uncertainties in 
the gas mass and distribution, we use S g = 10 g cm" 2 . The 
magnetic field strength in Arp 220 is not known, but it is 
very large. Estimates range from the minimum energy value 
of 0.3 m G to 30 mG, the maxi m um allowed by hydrostatic 
balance dThompson et al.1 120061) . iTorresi d2004l) found that 
B w mG in the main disk, but B«5 mG in the starburst nu- 
clei. We run a grid of models in B for values of 0.25 mG to 
16 mG for each nucleus. 

Unlike the Galactic Center, NGC 253, and M82, Arp 220 
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FIG. 7. — The CR electron (left) and photon (right) spectra of M82 for models with p = 2.2, fdjff(GeV) = 10 7 yr, and 7/ = 0.1. The line and shading styles are the 
same as Fig. [3] Synchrotron emission can dominate IC emission at energies less than a few keV in high B models. Data plo tted on rig ht are Fermi and VERITAS 
(triang les), as well as the dif fuse emission from Chandra (square) and total emission from BeppoSAX (Cappi et al. 1999), Suzaku iMiyawaki et al. 2009), and 
Swift (Cusumano et al. 2010) (stars). X-ray luminosities are scaled to one In bin in energy assuming Y = 2.0. 



has no 7-ray detection in either GeV or TeV bands. For a 
~ GeV upper limit, we use the > 100 MeV fl ux limit for 
5a sou r ces in the Fermi-L A^T one year catalog ( Ab do et alj 
l2010bl) . lAlbert et alj d2007l) provide upper limits in the VHE 
range for Arp 220 (specifically, an integrated photon flux from 
0.36 to 1.8 TeV). To keep the parameter space manageable, 
the CR proton acceleration efficiency 77 is assumed to be a 
standard value of 0.1. 

The X-ray emission o f Arp 220 has been studied by Chan- 
dra and XMM-Newton dlwasawa et alj|200ll IClements et alj 
12001 iMcDowell et al.1 120031: llwasawa etaill2005l) . XMM- 
Newton has detec ted a Fe K line, sugg esting at least some 
thermal emission (Iwasa wa et alJ 120051) : but there is also a 
hard X-ray continuum of unknown origin. There are sev- 
eral diffuse components, including a nuclear source roughly 
cospatial with the western nuclear starburst (Arp 220 X-l), a 
"hard halo" with a radius of a few hundred parsecs, a circum- 
nuclear halo extended on kiloparsec scales, an d vast "plumes" 
stretching out 10 kpc from th e galactic center (IClements et alj 
l2002t IMcDowell et al.ll2003h . The central unresolved source 
Arp 220 X-l is approximately coincident with the western 
nucleus, with an una bsorbed 2-10 keV l uminosity estimated 
at 4 x 10 40 ergs s" 1 (IClements et alJl2002l) . A second X-ray 
source, Arp 220 X-4 (L 2 -io » 1.5 X 10 40 ergs s" 1 , after cor- 
recting for absorption), is tentatively identified with the east- 
ern nucleus, although the geometry does not exactly match 
the radio morphology. For this paper, we assume that Arp 
220 X-l is the western starburst and Arp 220 X-4 is the east- 
ern starburst; in Table [5] we list the (unabsorbed) 2-10 keV 
synchrotron luminosities so they can be compared with other 
values for the luminosity. 

Modelled Radio and -f-rays - We show an example radio 
spectrum for the western and eastern nuclei in Figure [8] Sec- 
ondaries dominate the GHz synchrotron emission for B > 
4 mG, while primaries dominate for the lower magnetic field 
strengths. As with M82 and NGC 253, the radio data for each 
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FIG. 8. — Our predicted s ynchro t ron ra dio spectrum (red, solid) compared 
to observations compiled in Tones 1 2004). The open circles are for the west- 
ern nucleus and the open squares are for the eastern nucleus. The line and 
shading styles are the same as in Figuref2] 

nucleus appears to be flatter than predicted in our models. The 
discrepancy at ~ 10-20 GHz becomes larger in low B mod- 
els, since we are then probing the higher energy (more IC and 
synchrotron cooled) parts of the spectrum. However, it 
is possible that a thermal emission component flattens the ra- 
dio spectra. The thermal dust contribut ion is expected to be 
important only at 100 GHz and above dDownes & Solomon] 
ll99ot lTorres 2004). Furthermore, the 15 GHz data points ap- 
pear systematically high for b oth the west a nd e ast nuclei, and 
are high even in the models ofQbrjil d200j). As|Torres|([2004) 
note, the radio data use different beam sizes at different fre- 
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quencies. Finally, it is possible that bremsstrahlung and ion- 
ization are even stronger than expected. In order for advective 
escape to be strong enough to flatten the spectrum, wind es- 
cape times would have to be < 1000 yr (compare eqns.[5]and 
\12\i . requiring speeds of 3> 10 4 km s , which is unlikely. 

The non-detection of either nuclear starburst by Fermi al- 
ready implies that B > 250 pG in each of them, otherwise the 
leptonic emission alone would be observable. Note that we 
applied the GeV constraint to each nucleus independently; 
considering them together would push the lower limit on B 
towards ~ 500 pG. With more data, Fermi will be able to 
constrain the leptonic 7-ray emission further and increase the 
lower limit on B. Furthermore, note that in the B = 500 /iG 
and B = 1 mG models, we require that more power goes into 
CR electrons than protons. This is unlikely, considering that 
the reverse situation holds for the Milky Way. 

In all of our models, the nuclear starbursts of Arp 220 are 
proton calorimeters, with F ci \ = 0.84-0.98. Even at energies 
above 10 TeV, where diffusion is strongest, F c - d i « 0.6 when 
f diff (3 GeV) = 1 Myr. 

Synchrotron X-rays - In the extremely dense gas of Arp 
220, the CR protons at 10 - 100 TeV are efficiently converted 
into pionic products, including secondary e and 7-rays. In 
turn, the extreme FIR radiation field converts 10 - 100 TeV 7- 
rays into pair e ± . Finally, the secondary and tertiary e ± emit 
prodigious synchrotron X-rays. Models with p w 2.0 and slow 
diffusive escape can explain most of the observed X-ray emis- 
sion from the nuclei as synchrotron. Without any significant 
TeV 7-ray constraints on p, we cannot rule out these mod- 
els: it is possible that all of the hard X-ray emission from Arp 
220 is synchrotron (though, see the later caveats about X-ray 
absorption). 

For p = 2.2, the allowed synchrotron fractions are 0.9-42% 
(fiducial: 15%) for the western starburst and 2-84% (fiducial: 
34%) for the eastern starburst. In high B models (~ 4 mG) 
with p = 2.2, the synchrotron fraction is ~ 10% for the west- 
ern starburst and ~ 30% for the eastern starburst, consistent 
with previous expectations that HMXBs dominate the hard 
X-ray emission of starbursts. However, synchrotron is not so 
minor that it can be ignored, especially considering the uncer- 
tainties. Lower B models have a larger synchrotron fraction if 
7max = 10 9 and a smaller fraction if 7™™ = 10 6 . At still higher 
p = 2.4, synchrotron X-ray emission is reduced to a few per- 
cent of the observed diffuse X-ray flux in each starburst. 

We find that X-ray synchrotron emission in Arp 220's nu- 
clear starburst is harder than in the other starbursts, with 
T2-10 ~ 1.84-2.23. The hardest synchrotron spectra arise in 
low B models, where the Klein-Nishina bump is very promi- 
nent. With B greater than in M82 and NGC 253, lower elec- 
tron energies are being probed (grey shading in left panel of 
Figure|9]l, where the cutoff in the proton spectra does not mat- 
ter as much. To compare, the spectral slope of the diffuse 
hard X-ray emission in Arp 220 is not well constrained ob- 
servat ionally, but is p r obably harder than synchrotron emis- 
sion: | Clements et alj (120021) measured V = 1.4 ± 1, while 
llwasawa et al.l (120011) meas ured T = 1 .7 f o r the total X-ray 
emission of Arp 220, and llwasawa et al.l (120051) measured 
T = 1 .2+j} 4 , though for a Galactic absorption column. 

Inverse Compton X-rays - The contribution of Inverse 
Compton emission depends on the magnetic field strength. 
The contributions are highest in the low B models: for B = 
500 fiG, IC is 1.6 - 2.8 times the unabsorbed 2-10 keV lumi- 
nosities in the western starburst and 3.1 - 5.7 times the eastern 



starburst's luminosity. However, the actual X-ray absorption 
is uncertain, so these cannot be said to constrain the magnetic 
field with any certainty. At the other extreme of B = 4 mG, IC 
makes up ~ 6-7% (western) and 12- 14% (eastern) of the 
unabsorbed 2-10 keV luminosity. Thus, IC emission in our 
mod els can be much mo re inefficient than usually expected 
(c.f. llwasawa et al. 2001). This is because B is allowed to be 
much higher than the usual minimum-energy estimate; with 
the higher magnetic field strengths, IC is much weaker when 
scaling from the synchrotron radio emission. 

The synchrotron luminosity can be anywhere from 0.05% 
to 57 times the IC emission in Arp 220's western nucleus and 
0.06% to 35 times the IC emission in Arp 220's eastern nu- 
cleus, an even greater range than in the other starbursts be- 
cause p is not yet constrained by TeV 7-ray observations. 
Restricting ourselves to p = 2.2 models narrows the range 
(0.4-240% for the western nucleus; 0.6-280% for the east- 
ern nucleus), with fiducial values for the synchrotron-to-IC 
ratio of 2.3 for the western nucleus and 2.7 for the eastern 
nucleus. 

Like the other starbursts, IC emission is predicted to be 
harder than synchrotron emission, with T2-io ~ 1.30- 1.69, 
which is a closer match to the observed X-ray emission. It 
is possible, then, that Arp 220 has low magnetic fields of less 
than a milliGauss and its X-ray emission is IC emission. How- 
ever, both HMXBs and thermal emission are also expected to 
have the same spectral shape and dominate the X- ray spec- 
trum (iPersic & Rephaelill2002l: llwasawa et ai1l2005l) . 

Conclusion - Arp 220 is extremely efficient at converting 
VHE CR proton energy into synchrotron X-rays. Without 
the TeV 7-ray detections that exist for the Galactic Center, 
NGC 253, and M82, we cannot rule out hard spectra extend- 
ing to very high energies. Therefore, it is entirely possible to 
generate or even exceed the observed absorption-corrected X- 
ray luminosities of Arp 220 X-l and Arp 220 X-4 (assumed 
to be the western and eastern nucleus, respectively). There- 
fore, if the unabsorbed luminosities are correct, we can con- 
strain p = 2.0 electron spectra extending to multi-TeV ener- 
gies based on the X-rays alone. Even in p = 2.2 models, syn- 
chrotron can easily be ~ 10% or more of the X-ray luminosi- 
ties of each starburst. 

A large uncertainty in our models of Arp 220 is the amount 
of hydrogen absorption. A surface density of 10 g cm" 2 im- 
plies Nh ~ 6 x 10 24 cm" 2 , which is Compton thick. If this 
is correct, then the true X-ray luminosity of Arp 220 may be 
much higher. By contrast, the X-ray emission indicates hydro- 
gen column depths of Nh ~ 3 x 10 22 cm" 2 (T, g = 0.05 g cm" 2 ; 
IClements et alj [2002h . It is possible, however, that CR e ± 
emitting synchrotron X-rays are not cospatial with the nuclear 
starbursts themselves, either instead diffusing out of the nuclei 
(though this is unlikely g iven the extremely rapid losses in the 
extreme environments; IV5igfT989l:ICondon et al.lll991l) or are 
pair-produced from 7-rays outside the nuclei, in which case 
hydrogen absorption may be relatively small. It is also pos- 
sible that the covering fraction of the hydrogen is less than 
unity, in which case, some X-ray emission will be able to 
escape. Hydrogen absorption will also affect other explana- 
tions of the diffuse X-ray flux in Arp 220, such as IC emis- 
sion and high-mass X-ray binaries. Indeed, it might explain 
why Arp 220 seems underlum inous in hard X-rays for its star- 
formation ^ate J _asjiote(I_bjJlwasawaetay d2005l) and sup- 
ported bv lLehmer et all d2010h . 

5. DISCUSSION 



19 




FIG. 9. — The CR electron (left) and photon (right) spectra of Arp 220 for models with p = 2.2 and ?diff(3 GeV) = 10 7 yr. The line styles and colors are the 
same as in Figure \5\ The effects of neutral hydrogen absorption are not shown on right. Observed absorption-corrected luminosities from Chandra are plotted, 
scaled to one In bin in energy assuming T = 2.0. Synchrotron emission dominates IC emission at energies less than ~ 10 keV in high B models. 



5.1. Synchrotron, the FIR-X-ray Correlation, and 
Submillimeter Galaxies 

Star-forming galaxies show a correlation between their FIR 
luminosities and their 2-10 keV X-ray luminosities, with 
L2-10 ~ 1 0~ 4 Lfir when AGNs do not contribute to the X-ray 
emission (David et al. 1992; Grimm et al. 2003; Persic et al. 
12004 iPersic & Rephaelil 120071; iLehmer et all I2010I) . The 
correlation holds better for the most intense starbursts like 
ULIRGs, where there is less infrared and X-ray em ission 
from older stellar populations ( Persic & Rephaeli 2007). X- 
ray lumino sity has often been used as a star-formation rate 
indica tor dRanalli et a l. 2003; Grimm et al. 2003; Pers ic et ail 
120041) . The X-ray luminosity is often attributed to HMXBs, 
compact objects accreting matter from a massive companion. 
HMXBs are expected to trace the recent stellar population. Is 
it possible that synchrotron accounts for much of the X-ray 
luminosity in some galaxies? 

Our models show a trend, with increasing synchrotron X- 
ray emission with increasing gas density: synchrotron is a tiny 
fraction of the Galactic Center's X-ray luminosity if primaries 
cut off at low energies, a small portion of M82's diffuse hard 
X-ray luminosity, and possibly a substantial fraction of Arp 
220's hard X-ray luminosity. Three basic quantities affect the 
luminosity from secondary and pair e ± : (1) the amount of gas, 
which CR protons collide with to create pionic secondaries; 

(2) the amount of star-formation, which creates IR photons 
that pair produce off VHE 7-rays, and which is generally cor- 
related with the amount of gas through the Schmidt Law; and 

(3) the physical size of the starburst. 

At the extreme end of synchrotron X-ray luminosity from 
secondary and pair e^, it is worth considering the diffuse syn- 
chrotron emission from submillimeter galaxies (SMGs), well- 
studied extreme starbursts at high redshift (z > 2). SMGs con- 
tain large amounts of gas, so they should be relatively effi- 
cient at converting energy in CR protons into 7-rays and sec- 
ondary e ± . Furthermore, they are much more extended spa- 



tially than the more compact starbursts observed at low red- 
shift ( Chap man et al.ll2004t iTacconi et al.ll2006t iGenzel et al.l 
l2008tlLaw et al.ll2009H Younger et al.l l2010t) . This means that 
both diffusive and advective escape will be less effective at 
transporting CR protons out of SMGs, simply because they 
have to travel farther, further supportin g the case for pro- 
ton c alorimetry at multi-TeV energies dLacki & Thompson! 
2010). Finally, the radiation energy density in SMGs is very 
high, and because the SMGs are "puffy" with long sightlines 
through them, the optical depth to VHE 7-rays will be high. 
Therefore, we may expect that SMGs, like Arp 220, will be 
very efficient at converting energy in VHE CR protons into 
synchrotron X-ray emission with L^ nch > 10~ 6 Lsf- We can 
evaluate how good S MGs are at creating secondaries with the 
estimates in § 12.1.21 ITacconi et al.l d2006) find that SMGs have 
typical gas surface densities of Eg « 0.4 g cm" 2 . Using a typi- 
cal midplane-to-edge scale height h of 1 kpc, we find densities 
of (n) ~ 40 cm" 3 , comparable to those found by lTacconi et al.l 
(2006). From eqn. [TT] the typical pionic loss time in SMGs 
is fpion ~ 1-3 Myr. Eqn. [TJ] gives us advective loss times 
of fwind -3.0 Myr (/!/kpc)(v wind /300 km s" 1 )" 1 . The ratio 
of these times suggests that SMGs are proton calorimeters, 
depending on the unknown diffusive escape time and to a 
lesser extent on wind speeds (though they would have to be 
> 700 km s" 1 to reduce F ca \ significantly). 

To study this effect, we modeled generic galaxies of 
several gas surface densities from E s = 0.01 g cm" 2 to 
30 g cm" 2 . The magnetic fields were chosen to scale as B = 
6 /iG(E g /0.0025 g cm" 2 ) 07 , in line with our pre vious work 
on th e radio emission of star-forming galaxies (lLacki et al.l 
2010). The st ar-formation rate was calculated using the 
Schmidt law in iKennicutl (119981) . The interstellar radiation 
field was taken to consist of the CMB, a T = 10000 K grey- 
body field corresponding to the unobscured starlight from 
young stars, and a greybody corresponding to dust-obscured 
starlight (40 K in the compact starbursts and 30 K in the puffy 
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FIG. 10. — Fraction of the bolometric flux from star-formation Lsf in 2 - 
10 keV synchrotron (black/grey) and IC (blue) emission. On the bottom, we 
show h = 100 pc; the top shows h = 1000 pc, as is the case for SMGs. The 
lines indicate p = 2.2, £ = 0.01, and ?diff (3 GeV) = 10 Myr for compact star- 
bursts and f,jjff(3 GeV) = 100 Myr for puffy starbursts. Solid is 7^"™ = 10 6 , 



while dashed is 75 



10 . Shading indicates possible values for other val- 



ues of those parameters. For a typical star-formation hard X-ray luminosity 
Z.2-10 = 10 Lsf, synchrotron and IC are subdominant; for optimistic choices 
of parameters, synchrotron can be 10~ 5 Lsf. The circles represent the diffuse 
X-ray emission for the studied starbursts, while the stars represent the total 
(point so urces includ ed) X-ray emission that would be seen at large distances 
(from the lCappi eta l. 1999 fluxes for M82 and NGC 253, scaled to our dis- 
tances; from lMcDowell et aU20"03l for Arp 220). 

starbursts). Scale heights of h = 100 pc and 1000 pc were 
modeled. We chose r\ = 0.1 and varied £ from 0.003 to 0.03; 
several values of fdiff(3 GeV) and 7m™ were tried. 

Our models give the ratio of both synchrotron and IC emis- 
sion to the bolometric power from star-formation at 2-10 
keV, and these are plotted in Figure [10] Both components 
are always a minority of the standard hard X-ray luminos- 
ity, 10~ 4 Lsf! from the FIR-X-ray correlation. Both the syn- 
chrotron and IC fraction increase with E g . IC shows less 
variation in its flux, since it mostly comes from the relatively 
well-constrained low energy CR electron spectrum, although 
L2-10 still varies by a factor of 3-35 for a given T g when 
h= 100 pc (5-180 for h = 1000 pc). At the highest E g , the IC 
emission is 2% or less of the typical 2-10 keV luminosity for 
h = 100 pc (6% or less when h = 1000 pc), falling to < 0.1% 
for T g = 0.01 g cm" 2 (blue shading). 

Synchrotron X-ray emission shows a much greater vari- 
ation with different parameters (see grey shading), ranging 
from < 0.1% (assuming p = 2.4 or that T g is small) to 24% of 
the 2-10 keV luminosity of typical starbursts. Models with 
high are much more efficient synchrotron radiators at 
low Eg. At high Tig, however, the synchrotron fraction starts 
to converge, as CR proton energy is efficiently converted to 
pions and 7-ray photons are efficiently converted to e^- pairs. 
For our fiducial values of p = 2.2 and f m ff(3 GeV) =10 Myr, 
^40* approaches 4% of the total 2-10 keV luminosity in the 
densest starbursts when h = 100 pc (3% when h = 1000 pc). 
Synchrotron dominates IC for our standard parameters when 
^ 0.1 g cm" 2 and h = 100 pc, which are typical of com- 
pact starbursts. Thus, we do not expect synchrotron emission 



to usually cause large deviations in the X-ray luminosity, al- 
though it is present, and at greater levels than IC. 

As with the more compact starburst models, the syn- 
chrotron fraction of 2-10 keV emission of the puffy starburst 
models varies over many orders of magnitude with param- 
eters. For our fiducial f<jiff(3 GeV) = 10 Myr from compact 
starbursts, we actually find that the synchrotron fraction is de- 
creased in puffy starbursts when 7 m '™ is small. This is be- 
cause the diffusive escape time is the same, while the density 
is 10 times lower at a given T g ; the lower pionic yield reduces 
the amount of secondary e ± and the pionic 7-rays that seed 
pair e^. However, for a given diffusion constant, diffusion out 
of the larger SMGs takes longer. Setting f diff (3 GeV) to 100 
Myr restores the synchrotron fraction to near its compact star- 
burst value. Comparing models with fdiff(3 GeV) =100 Myr 
for h = 1000 pc and f diff (3 GeV) = 10 Myr for h = 100 pc, and 

7111™ = 10 6 , there is a ~ 50% enhancement in the synchrotron 
emission for E g = 0.03-0.3 g cm" 2 , because the 77 optical 
depth is greater along the long sightlines in puffy starbursts. 
Given the uncertainties in the parameters and the small contri- 
bution of synchrotron at these energies, this is not a big effect. 

An interesting effect that occurs in our h = 1000 pc mod- 
els is that the IC fraction actually increases. Unlike the e ± at 
multi-TeV energies, the e ± responsible for 2-10 keV IC emis- 
sion are relatively low energy (~ 100 MeV- 1 GeV, emitting 
synchrotron at ~ GHz frequencies). In compact starbursts, 
these e ± are cooled significantly by bremsstrahlung and ion- 
ization losses, suppressing their IC emission. In puffy star- 
bursts with the same T g , however, the gas volume density is 
lower, so that bremsstrahlung and ionization losses are not as 
important. This means that there is more power to go into 
IC (and possibly synchrotron) losses at these energies, en- 
hancing the X-ray IC (and possibly GHz synchrotron) emis- 
sion in puffy starbursts rel ative to compact starbursts (see also 
lLacki & T hompson 201(J. 

In fact, intense starbursts like SMGs, have relatively weak 
X-ray emission compared to other starbursts, if t hey lack an 
AGN, with L x « IO^Lsf (lAlexander etalJ l2005). This sug- 
gests that diffuse synchrotron accounts for a few percent of 
the total hard X-ray emission from SMGs, and not just the 
diffuse hard X-ray emission. 

5.2. What Neutrinos and TeV 7-rays Can Tell Us 

There are many uncertainties in our models, leading to a 
broad range in the synchrotron X-ray luminosities. These in- 
clude the hardness of the e ± spectrum, the strength of B which 
sets the importance of secondary e ± , and the efficiency of 
pionic losses at multi-TeV energies. The synchrotron X-ray 
emission is enhanced significantly in starbursts where pionic 
secondary e ± and pair production e ± are produced at a high 
rate, especially if p ss 2.0. A necessary byproduct of pionic 
e and 7-rays is pionic neutrinos. Starbursts that are bright in 
neutrinos should therefore be relatively bright in synchrotron 
X-rays. Of course, the converse is not necessarily true if 
the synchrotron X-ray luminosity comes largely from primary 
e . TeV 7-rays are also generated by both high energy pro- 
tons and electrons in the intense environments of starbursts. 
While ~ 10-100 TeV 7-rays are heavily absorbed by pair- 
production, there is a window at lower energies below a few 
TeV (see Figure[T]i. Indeed, TeV 7-rays are the most powerful 
constraint yet on synchrotron X-ray emission from NGC 253 
and M82. 

M82 is in the Northern hemisphere and visible to IceCube. 
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FIG. 1 1 . — The all-flavor neutrino (plus antineutrino) spectrum of M82 in 
various models with r\ = 0.10. Grey is p = 2.0 and black is p = 2.2; solid is 
*diff(3 GeV) = 1 Myr, dashed is 10 Myr, and dotted is 100 Myr. The neutrino 
emission directly scales with rj. 

In Figure QT| we see that the predicted neutrino spectrum of 
M82 varies widely between different models. Models with 
low p and high t^s naturally predict that M82 is very bright in 
VHE neutrinos: the more protons at high energies, the more 
neutrinos are produced. The magnetic field strength B is also 
linked to the neutrino flux through the radio and 7-ray ob- 
servations. In models with high B, the observed radio flux 
is mostly secondary e^, and the accompanying pionic 7 rays 
make up almost all of the detected GeV to TeV 7 rays (as seen 
in Figure |7). The 7-ray flux therefore translates into a high 
neutrino flux. However, in models with low B, the observed 
radio flux is mostly primary e^, and leptonic emission makes 
up a majority of the GeV (and possibly TeV) 7-ray flux. In 
order not to overproduce the 7-rays, the CR protons must be 
less efficiently accelerated, in turn implying lower neutrino 
flux. 

The muon neutrino flux of M82 at 10 TeV is anywhere be- 
tween 1(T 14 and 10" 12 TeV cm" 2 s" 1 for p = 2.0-2.2. Un- 
fortunately, M82 is at high declination (sin S sa 0.93), wher e 
IceCube's sensitivity is relatively weak ( Abb asi et al.l l2009a). 
Even with a full year of IceCube-80, our most optimistic mod- 
els predict that M82 is ~ 5 times too faint to be observed. 
NGC 253 would be even fainter, and since it is located in the 
Southern Hemisphere sky (sin 5 ~ -0.43), the available ex- 
isting neutrino dete c tors are far less se nsitive (lAbbasi et al.1 
2009b; ICoyie1l20Tl lAbbasi etalJl201lh . However, the dif- 
fuse neutrino background from starbursts may be detectable 
by IceCube, w hich may shed some ligh t on the CR spectrum 
at high energy (|Loeb & W axman 2006). Stacking searches of 
nearby, bright star bursts may also lead to a detectable signal 
dLacki et alJl20TTb . 

Arp 220 is much fainter from Earth than M82, far beyond 
the reach of neutrino detectors like IceCube, but currently 
has no TeV 7-ray detection (th e current upper limit is set by 
MAGIC in I Albert et all 120071) . Thus, high p is allowed by 
the data. Just as TeV 7-rays constrained M82, TeV 7-rays 
can constrain the fraction of Arp 220's X-ray emission that 
comes from synchrotron. In FigureQ~2] we show the predicted 
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FIG. 12. — The Earth-observed pair production absorbed VHE 7-ray spec- 
trum of Arp 220 in various models, for 7ma 1 x 

= 10 6 (black for p = 2.2 and grey 
for p = 2.0) and 7^"™ = 10" (dark blu e for p = 2.2 and light blue for p = 2.0). 
The line styles are the same as in Fig. 1111 

VHE 7-ray flux of Arp 220's western and eastern nuclear star- 
bursts in several models compared to the expected 5a sensitiv- 
ity of VERITAS and Cheren kov Telescop e Array (CTA) after 
50 hours of integration time (Doro 2009|). 

VERITAS will not be able to detect Arp 220 at 5 a ex- 
cept in our most optimistic models with high 7m™, low B, 
and p = 2.0 (Figure [12] lower panel). However, if B is low, 
the proposed CTA will be able to detect Arp 220 except in 
more pessimistic models with p = 2.2 and small 7m™. If B 
is high, even CTA will have trouble detecting Arp 220's nu- 
clear starbursts (FigureQ~2j upper panel), especially if p = 2.2. 
We note, though, that our fiducial supernova rate for both of 
these starbursts (0.7 yr -1 ) is about one-third of that expected if 
the entire TIR luminosity of Arp 220 is due to star-formation 
(expected supernova rate of 2 yr" 1 ). Thus, the total TeV lu- 
minosity of Arp 220 may be ~ 3 times larger, in which case 
CTA will be able to detect it even in high B cases. Both the 
synchrotron X-ray flux and the TeV 7-ray flux are insensitive 
to the diffusive escape time: the protons are being efficiently 
trapped and creating secondaries at TeV energies. The effects 
of pair production absorption are clearly visible in Figure [121 
the 7-ray flux should plummet at energ ies higher tha n a TeV, 
an effect first predicted for Arp 220 by lTorres I(l2004l) . 

6. CONCLUSION 

Starbursts are luminous in hard X-rays and accelerate large 
amounts of CRs. We have explored whether CR e in star- 
bursts could generate the observed diffuse hard X-ray emis- 
sion through synchrotron emission. The strength of this emis- 
sion depends on the very poorly constrained 10 - 100 TeV CR 
e ± spectrum in starburst galaxies. Using one-zone models to 
predict this emission, we also reconsider the previously sug- 
gested contribution from IC emission of 0.1 - 1 GeV CR e , 
in light of the magnetic field strengths expected in starburst 
galaxies. Our conclusions are as follows: 

• We have considered the energetics for synchrotron X- 
ray emission in section [2] with simple order of mag- 
nitude estimates. The diffuse hard X-ray emission of 
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M82 could be synchrotron if CR escape is not important 
at PeV energies and the CR spectrum is hard enough 
(p s» 2.0). Much of this emission would come from the 
previously neglected pair-production e ± , which are ef- 
ficiently generated from pionic 7-rays in the intense IR 
bright environments of starbursts. 

• We have also constructed one-zone models of several 
starbursts with standard spectra to compare with obser- 
vations in radio and 7-rays. The magnetic field strength 
in these models is limited at the low end by constraints 
on leptonic 7-ray emission and at the high end by con- 
straints on radio emission from secondary e , For 
the physical conditions we assume (listed in Table [TJ, 
these limits are 50 fiG < B < 150 ^G in NGC 253, 
50 < B < 200 [iG in M82, and 0.5 mG < B < 4 mG 
for Arp 220's nuclei. In low B models there is a large 
component of primary electrons, while in high B mod- 
els the secondary <? ± dominate at GeV energies in M82 
and Arp 220's nuclei. 

• The existing 7-ray data of the modelled starbursts rules 
out most p = 2.0 models, and we find that the syn- 
chrotron emission is probably a minority of the hard 
X-ray emission in the Galactic Center, NGC 253, and 
M82 with our assumptions. Nonetheless, there is still 
enough variation in the parameters for synchrotron to 
make up from less than one to several tens of percent 
of the synchrotron emission. The synchrotron X-ray 
emission is highest in low B models with large 7^™, 
lower in high B models, and lowest in low B models 
with small 7™™. In our fiducial models (with high 

B and large 7m 1 ™), the fraction of unresolved hard X- 
ray emission contributed by synchrotron is 6% in the 
Galactic Center, 9% in NGC 253, 2% in M82, 15% in 
Arp 220's western nucleus (Arp 220 X-l), and 34% in 
Arp 220's eastern nucleus (Arp 220 X-4). These re- 
sults are consistent with previous predictions that X-ray 
binaries contribute most of th e X-ray emission in star- 
burst galaxies at a few keV dPersic & Rephaelil 120021 : 
iPersic et al.ll2004h . We find that the 2 - 10 keV X-ray 
synchrotron emission is 10~ 7 - 10~ 6 Ltr in these fiducial 
models; in models of generic star-forming galaxies it 
peaks at - 4 x 10" 6 L IR for p = 2.2. 12 

• Pair production is predicted to contribute significantly 
to the high energy e ± population in NGC 253, M82, and 
Arp 220, with pair e^- comparable to the pionic e ± . The 
huge densities of IR starlight photons efficiently con- 
vert ~ 10— 100 TeV 7-rays into e ± pairs. These pairs 
then can efficiently radiate X-ray synchrotron. How- 
ever, the actual density of pair e is directly related to 
7-ray emissivity at very high energies, which in turn de- 
pends strongly on the number and escape times of CR 
protons at these energies. The Galactic Center region is 
transparent to VHE 7-rays, so that pair e are negligi- 
ble. 

• We find that Inverse Compton emission in the 2-10 
keV band is also a minority of the X-ray emission. In 
our models, the magnetic field energy density can be 

12 For comparison, according to the observed FIR-radio correlation, the 
GHz radio luminosity is ~ l(T 6 Z,i R )Yun et aU200Tl) . 



greater than some of the values derived with the mini- 
mum energy estimate. This means the IC emission ex- 
pected by scaling the GHz radio emission is lower than 
in previous estimates. In our fiducial models, the frac- 
tion of unresolved X-ray emission contributed by IC is 
0.1% in the Galactic Center, 3% in NGC 253, 1% in 
M82, and ~ 10% in Arp 220's nuclei. The IC emis- 
sion at 2 - 10 keV is much harder (IVio = 1.2—1.7) 
than the synchrotron emission at the same energies 
(IVio = 1.8-2.5). 

• Key uncertainties include the maximum energy of pri- 
mary electrons (7™™), the rate of escape of CR protons 
at high energies (fdiff), and the spectral slope of CRs 
(p). However, some of these quantities might be con- 
strained by future neutrino telescopes, which can de- 
termine the pionic luminosities of starbursts at TeV to 
PeV energies. Future TeV 7-ray telescopes like CTA 
can constrain the electron spectrum at high energies as 
well, especially in Arp 220 which has relatively weak 
TeV upper limits. The effects of hydrogen absorption in 
dense starbursts also need to be explored, although this 
applies to IC emission and other sources of the X-ray 
flux, not just synchrotron. 

The large X-ray foregrounds from HMXBs and other 
sources makes detection of the synchrotron X-rays difficult. 
However, if the synchrotron emission can be detected, it will 
have important implications for our understanding of CRs in 
starbursts. Very little is known about the CR population at 
very high energies in starbursts, such as what the escape rate 
is or whether there is a "knee" in the CR spectrum as there is 
in the Milky Way. Combined with upcoming neutrino mea- 
surements, detection of synchrotron X-rays can extend our 
understanding from the directly observed TeV energies. 

Even in the Milky Way, we know relatively little 
about the CR electron population at multi-TeV energies 
dKistler & Yuksell |2009), and the propagation of electrons 
at 100 GeV - TeV energies is still poorly understood 
(IChang et all 120081: lAdriani et alJ 120091 lAbdo et alJ 120091 
Aharonian et al. 2009a). It is possible that we are underes- 
timating the primary population significantly. We assumed 
that the primary CR electron injection spectrum is a power 
law extending to TeV-PeV energies, with one source such as 
supernova remnants dominating the primary electron spec- 
trum. An additional hard primary component (from PWNe, 
for example) on top of a softer component that explains the ra- 
dio synchrotron emission would increase the synchrotron X- 
ray emission. However, a second component such as PWNe 
would have to be very hard (p < 2.0) to enhance the X-ray 
flux significantly; otherwise IC emission from the low energy 
e^ may exceed the observed 7-rays. 

If a significant fraction of the hard X-ray emission from 
starbursts is synchrotron, it may be polarized. Although the 
radio synchrotron emission of M82's halo is strongly po- 
larized, the starbu rst core itself has very little polarization 
dReuter et al.l 1994b . This may be because of Faraday depolar- 
ization, in which the back of the starburst is strongly Faraday 
rotated with respect to the front. Faraday depolarization, how- 
ever, will be completely unimportant at X-ray frequencies. 
The low polarization could also arise if the magnetic field was 
turbulent and isotropic on small scales, in which case the X- 
ray emission will also be unpolarized. However, even turbu- 
lent magnetic fields can be made anisotropic through shear 
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and compression (e.g., iLaind Il980t ISokoloffet all 119981) . 
Furthermore, there is observational evidence of anisotropic 
magnetic fields in M82 fr om submillimeter and infrared po - 
larization measurements (|Greaves et al.1 120001: 1 Jones! 2000). 
ISiebenmorgen & Efstathioul (120011) and feeiffert et all (1200% 
found that Arp 220 has a low infrared and submillimeter po- 
larization, though their source apertures blend the two nu- 
clear starbursts and the outlying starburst disk together. In 
any case, the diffuse hard X-r ay emission of M82 only has 
a flux of ~ 0.1 -0 2 milliCrab dStrickland & Heckmanll2007t 
iKirsch et al.ll2005l) : only if it had large polarization could the 
Gravity and Extreme Magnetism Small Explorer (GEMS) de- 
tect its polarization (Swank et al. 2009). 

It is tempting to consider even higher energy synchrotron 
emission, in the 10 keV - MeV band, to directly probe 30 TeV 
- PeV e ± . Presently, there are only relatively weak upper lim- 
its (and a claimed detection of NGC 253) of this emission of 
nearby starbursts and ULIRGs from the OSSE instrument on 
the Compton Gamma-Ra y Observatory (Bhattachary a et ail 
Il994t iDermer et alll 19971). and a re cent detection of M82 by 
Swift-BAT dCusumano et al.l [2010). With its much greater 
sensitivity in the 10-80 keV range over previous tele- 
scopes, NuS TAR can improve this situation (lHarrison et al.1 
120051 120101) . However, above ~ 10 keV, we generally ex- 
pect the very hard IC emission to finally bury the softer syn- 
chrotron emission, even if the CR electron spectrum does ex- 
tend to multi-PeV energies (compare the synchrotron and IC 
emission in M82 in Figure [7] right panel). Furthermore, the 
angular resolution of NuSTAR is 40", corresponding to 750 pc 
at the distan ce of M82, so point sourc es may contaminate 
the emission (lHarrison et al.ll2005l 12010). On the other hand, 
even a measurement of the IC emission itself could be useful 
in constraining the magnetic field strength. The synchrotron 
luminosity may remain within a factor of ~ 2 of the IC until 
MeV energies in strong starbursts like Arp 220, if they have 

13 Such MeV emission would be unaffected by hydrogen absorption as 
well, which would be ideal for an extremely dense starburst like Arp 220, 



very strong magnetic fields (B > mG; as seen in Figure [9j. 13 
Even if the IC and point source emission was not a signifi- 
cant foreground, ultimately the synchrotron will be buried by 
bremsstrahlung and finally direct pionic 7-ray emission above 
10 MeV. Since there is more power in lower energy protons 
than VHE protons, the synchrotron emission cannot exceed 
the pionic emission at higher energies. Thus, synchrotron can- 
not be used to probe the CR e ± spectrum above about a PeV, 
unless an enormous new component of primaries is present. 

While we have considered mainly 7-rays from star- 
formation, 7-rays in starburst environments can come from 
other sources: most notably, an AGN. In particular, Sgr A* 
in o ur own Galaxy is kno wn to be a source of VHE 7- 
rays dAharonian et al.ll2004l) . It resides in the Central Clus- 
ter, a dense star cluster with a FIR energy density of ~ 
6000 e V cm -3 over a region ~ 2 pc i n diameter (fTblesco et all 
Il996t see also lDavidson et alJIl^lHopkins et alj|2010h . In- 
terestingly, the VHE emission appears to hav e an exponen- 
tial cutoff at 15 TeV dAharonian et alj|2009bh . similar to the 
Klein-Nishina threshold for pair production on the FIR light 
of a starburst. Diffuse hard X-rays are observed from the re- 
gion near Sgr A*; these potentially could have a synchrotron 
contribution from pair e . Similar considerations may apply 
to other central stellar clusters around VHE-emitting AGNs. 
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although it is beyond the energy range of NuSTAR. 



APPENDIX 

OTHER PIONIC CROSS SECTION AND LIFETIME PARAMETRIZATIONS 

The spectrum of pionic secondary 7-ra ys and e at TeV energies depends sensitively on the behavior of the pion production 
process at these energies dKarlssonll20 08). The physics of pion production enters two ways in our models: (1) the lifetime of CR 
protons to pionic processes, which determines the steady-state CR spectrum and (2) the differential cross sections for production 
of pionic secondaries for protons of each energy. 

We run models where we vary the pionic lifetimes or the differential cross sections to consider the effects. We use the fiducial 
parameters listed in TableQ]for the Galactic Center, NGC 253, M82, and the two nuclei of Arp 220. We find the thermal fraction 
and £ through radio spectrum chi-square fitting (or radio flux normalization for Arp 220) described in § |3.3| for each variation on 
the pionic physics. We also compare the purely hadronic flux using the different assumptions of the pionic physics. 

Pion ic lifetime variations - In most of our models, we calculate the pionic loss lifetime by integrating up the iKamae et al.l 
d2006l) cross-sections for 7-ray, n eutrino, and e ± production (equation |33Tl. A commonly used p ionic loss lifetime comes from 
Mannheim & Schlickeiserl d 19941) (MS94) and is given in equationQT] Finally, Schlickeiser ( 2002) recommends a different pionic 
loss lifetime, which is longer at low energies but substantially shorter at VHE energies: 

t n = 2.2 x 10V 7" a28 (n/cirr 3 r 1 (Al) 

The new 7 dependence comes from the multiplicity of pions produced per collision. Formally, this equation is only valid for 
7 < 3 000, or below a few TeV. 

The Schlickeiser (2002) pionic lifetime predicts substantially less hadronic synchrotron X-ray emission, unless escape domi- 
nates the multi-TeV proton lifetimes. In our Arp 220 models, where pionic losses dominate even at these energies, the hadronic 
2-10 keV synchrotron flux is only 1/4 of the predictions using the cross-section derived t^. This is because the 7" 28 in 
becomes very small for 10 - 100 TeV CR protons, reducing the steady-state number of CR protons that are the source of sec- 
ondary and pair e ± . The secondary source functions do not correspondingly increase, so the pionic luminosity is underpredicted. 
The M82 model also gives only 64% of the hadronic synchrotron flux predicted with the cross-section integrated t^. The TeV 
7-ray flux is also reduced in these models by similar amounts. However, the synchrotron flux in the fiducial models of NGC 
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253 and the Galactic Center are affected by less than 10%, because diffusion sets the total number of protons in these models. 
Furthermore, the total synchrotron X-ray flux is barely affected in the Galactic Center, NGC 253, and M82, and only reduced to 
40% of nominal in the Arp 220 nuclei models, because primary CR e ± contribute heavily to the synchrotron emission. 

If we use the MS94 t^, the hadronic synchrotron X-ray emission is instead slightly enhanced. In Arp 220 models, the syn- 
chrotron X-ray emission is 113% that using our standard pionic lifetime. The total synchrotron X-ray emission is significantly 
greater in MS94 models, because higher £ are picked by our fitting processes: the synchrotron is 1.3 times greater in Arp 220 
models using MS94 and 1.2 times greater in M82 models. There are negligible differences in the TeV 7-ray fluxes between the 
MS94 models and eqn. 1331 models, but the GeV fluxes are only 80% as big in the M82 models and 70% as big in the Arp 220 
models. 

The[Schlickeiser (2002) pionic lifetime is only meant to be valid below 2.8 TeV, which is below the energy of the CR protons 
responsible for hard X-ray emitting e . Furt hermore, it i s deriv ed by integrating the energy of the pionic secondaries, much like 
we do to derive eqn. [33] Since we used the Kam ae et al.1 ([2006) cross sections in eqn. [33] which are explicitly valid up to several 
hundred TeV, we believe that eqn. [33] is more likely to be accurate in deriving the X-ray synchrotron emission. Thus, the pionic 
lifetime would only contribute less than a factor of ~ 2 uncertainty, much less than the uncertainties from 7S™, p, and f<jiff. 

Differential cross sec tion variations - There are several possible p arametrizations of the pionic cross sections. In addition to 
thelKamae et al.1 (120061) cross sections, th e GALPROP cro ss sections dMoskalenko & Stronglll998:IStrong & Moskalenkoll 19981: 
Stron g"et al.ll2000t) based on the work of iDermerl (Il986a|) are commonly used (in turn based on Steckerlll970l iBadhwar et al.l 
19771 and lStephens & Badhwarlll98U see also iDermerll 1 986bl) . 

Both the Ka mae et al.l (120061) and GALPROP cross sections are based on low er energy data, and IKamae et al.l (|2006) is only 
formally valid for CR protons up to 512 TeV. The pionic cross sec tions given in iKelner et ail d2006l) are valid for CR protons at 
energies at 100 GeV to 100 PeV. We consider a variation where the lKamae et al.l (120 06) cross sections are used for CR protons at 
energies below 100 GeV, and the lKelner et al.l (120061) cross sections are used above 100 GeV. 

Using the GALPROP cross sections reduces both the hadronic and total synchrotron X-ray luminosities. For the Galactic 
Center, the hadronic (total) GALPR OP synchrotron X-r ay luminosities are 67% (100%) the lKamae et al . (2006) luminosities; for 
NGC 253, they are 73% (83%) thelKamae et al.l (120061) luminosities; for M82, the GALPROP synchrotro n X-ray luminosities 
are 79% (75%) the iKamae etafl (120061) luminosities; and for Arp 220's nuclei, they are 105% (90%) the iKamae etail (120061) 
luminosities. Smaller £ are preferred using the GALPROP cross sections, reducing the leptonic contribution to the synchrotron 
X-ray emission. In the Arp 220 models, using the GALPROP cross sections leads to small increases of the hadronic 7-ray 
luminosities, about ~ 15% at GeV and ~ 2 0% at TeV energies. 

If we instead use the Kelner et al. ( 2006) cross sections for E p > 100 GeV, we find even less of an effect. For the Galactic Cen- 
ter, the hadronic (total) lKelner et all d2006l) synchrotron X-r ay luminosities are 84% (100%) the Kamae et al. (2006) luminosities; 
for NGC 253, they are 87% (98%) the lKamae etail ([2006) lumi nosities; for M82, th ey are 92% (97%) the lKamae etail (120061) 
luminosities; and for Arp 220's nuclei, they are 1 14% (110%) the Kamae et al. (2006) luminosities. The TeV 7-ray luminosity is 
27% brighter in the Arp 220 models using these cross sections instead of our standard cross sections. 

Thus, using other parametrizations of the pionic cross sections does not seem to add even a factor of ~ 2 uncertainty. 

OTHER MAXIMUM ENERGY PROTON CUTOFFS 

Throughout this work, we have assumed that primary CR protons are accelerated to a maximum Lorentz factor 7 max = 10 6 , or 
a maximum energy of ~ 938 TeV. This is loosely based on the observed knee in the CR spectrum in the Galaxy, but it is possible 
given current observational data that the knee is a t a different energy in starburst regions. Since the average energy of pionic 
electrons is about 1/20 of the proton energy (e.g.. IKelner et alJl2006l) . this translates to a pionic e ± cutoff near 50 TeV. Thus, 
models with low 7m™ or high B can be affected by the proton energy cutoff. To consider this effect, we rerun our fiducial models 
with 7 max of 10 s , 10 6 , and 10 7 . We consider both the total and the hadronic fluxes, as in Appendix lAl As in o ur fiducial models, 
we use the Ka mae et al.l (2006) cross sections for proton energies below 500 TeV and the lKelner et alJ (120061) cross sections for 
proton energies above 500 TeV. 

In models where 7 max = 10 5 , the hadronic synchrotron X-ray flux drops by almost an order of magnitude in most of the starbursts 
we consider. The hadronic synchrotron X-ray flux is only 9.4%, 8.5%, and 1 1% of its fiducial value in the Galactic Center, NGC 
253, and M82 respectively. In Arp 220's nuclei, the effect is not so severe since B is ~ 40 times higher than in the other starbursts, 
and synchrotron X-rays probe e ± energies that are ~ 6 times lower (eqn.Q]). However, the low energy proton cutoff still reduces 
the hadronic synchrotron X-ray flux to 43% of its fiducial value. The 2-10 keV hadronic synchrotron X-ray flux is also much 
softer in these models, because of the cutoff in the secondary e ± spectrum: the hadronic TJ^q i s 3.45, 3.29, 3.06, and 2.46 in 
the Galactic Center, NGC 253, M82, and Arp 220's nuclei. The total synchrotron X-ray flux is affected less than the hadronic 
flux, but note that we use 7m™ = 10 9 in our fiducial models so that primaries dominate in most of the starbursts. Thus the total 
synchrotron X-ray flux is 99.7%, 90%, 64%, 58%, and 49% of the fiducial values in the Galactic Center, NGC 253, M82, Arp 
220 West, and Arp 220 East, respectively. 

In models where 7 max = 10 7 , the hadronic synchrotron X-ray flux in the 2-10 keV band is enhanced by ~ 40%. The hadronic 
synchrotron X-ray flux is 141%, 141%, 136%, and 112% of its fiducial value in the Galactic Center, NGC 253, M82, and Arp 
220's nuclei respectively. The hadronic synchrotron X-ray spectrum is also hardened slightly: the hadronic is 2.32, 2.25, 

2.21, and 2.11 in the Galactic Center, NGC 253, M82, and Arp 220's nuclei. The total synchrotron X-ray flux is affected very 
little, with total synchrotron fluxes 100.2%, 103%, 114%, 110%, and 112% of fiducial in the Galactic Center, NGC 253, M82, 
Arp 220 West, and Arp 220 East, respectively. 

We conclude that if the secondary and pair e dominate at 10 - 100 TeV energies, the maximum proton energy can significantly 
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affect the synchrotron X-ray flux and spectral slope, particularly if the cutoff is much lower than a PeV. 
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Galactic Center X-ray luminosities 
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a Fraction of the observed diffuse 2-10 keV X-ray luminosity (scaled from Koyama et al. (1996) as described in § |4. l\ from that synchrotron emission accounts for. 

Ratio of the 2-10 keV synchrotron and IC luminosities. 
C Ratio of the predicted GeV 7-ray emission with the upper limit from EGRET in Hunter et al. 1 1997). 

^ Ratio of the predicted TeV 7-ray emission with the observed TeV flux observed from HESS in lAharonian et all 120061) and lCrocker et all 1201 lbl) . 

e This model satisfies our constraints (§ |3.3| . While r\ is not the minimum or maximum allowed or the closest to reproducing the TeV luminosity for one of the 7^""\ this is true for 

the other 7^"™ so we list it anyway, 
f 

Chosen as our fiducial model. 
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TABLE 3 
NGC 253 X-RAY LUMINOSITIES 
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2.37 


1.3 


0.66 




0.57 


0.081 


1.4e37 


0.016 


0.24 


2.37 


1.6 


1.8 


oo 


0.050 


0.15 


9.0e36 


0.011 


0.10 


2.33 


1.0 


0.53 




0.14 


0.13 


2.6e37 


0.030 


0.32 


2.33 


1.1 


1.5 


10 


0.0088 


















0.013 


















0.10 


0.14 


4.2e36 


0.0050 


0.051 


2.38 


1.1 


0.63 




0.28 


0.12 


1.2e37 


0.014 


0.15 


2.38 


1.4 


1.8 


1 


0.10 


















0.28 


















0.40 


0.12 


2.2e36 


0.0026 


0.027 


2.40 


1.4 


0.67 




0.57 


0.1 1 


3.1e36 


0.0036 


0.038 


2.40 


1.7 


0.94 


oo 


0.0088 


















0.050 


















0.20 


0.34 


4.8e36 


0.0056 


0.037 


2.36 


1.2 


0.59 




0.57 


0.23 


1.3e37 


0.016 


0.13 


2.36 


1.9 


1.7 


10 


0.0088 


















0.050 


















0.28 


0.31 


1.6e36 


0.0019 


0.013 


2.42 


1.3 


0.54 




0.57 


0.24 


3.2e36 


0.0038 


0.030 


2.42 


1.9 


1 . 1 


1 


0.0088 


















O.OjO 


















U.j / 
















J u 


n 1 A 


U.UZo 


A loll 
4. JCJ / 


ft CiKft 
U.UJU 


z.z 


9 11 
Z.J J 


i) ^i 
U.J J 


1 .o 


I 


n 1 a 
U. 14 


















U.zU 


u.uzy 


/.oejo 


u.uuyz 


U. Jo 


z. JJ 


U.J J 


U.oz 




0.57 


0.017 


2.2e37 


0.026 


1 . 1 


2.35 


0.81 


1 .7 


oo 


0.14 


0.035 


4.0e37 


0.047 


1 .0 


2.31 


0.54 


1 .4 




U.zU 


U.UJU 


^ loll 

j. /ej / 


U.Uo / 


1 1 
Z.J 


i i 1 
Z.J 1 


U.OD 


Z.U 


1 u 


n 1 A 
U. 14 


u.ujo 


y.oejo 


Ann 
U.U1Z 


n iq 
U.Jo 


o in 

Z.J / 


U.jZ 


U.oJ 




U.Zo 


U.UZ4 


z.uej / 


U.UZJ 


U.ol 


n in 
Z.J / 


u. /y 


1 . / 


1 


U.ZU 


















0.40 


0.025 


3.6e36 


0.0043 


0.14 


2.39 


0.85 


0.65 




0.57 


0.014 


5.2e36 


0.0061 


0.21 


2.39 


1.1 


0.92 


oo 


0.20 


0.068 


7.6e36 


0.0090 


0.23 


2.34 


0.62 


0.56 




0.40 


0.035 


1.5e37 


0.018 


0.51 


2.34 


1.0 


1.1 




0.57 


0.0064 


2.2e37 


0.025 


0.80 


2.34 


1.4 


1.6 


10 


0.28 


0.057 


2.7e36 


0.0032 


0.085 


2.41 


0.78 


0.53 




0.40 


0.039 


3.9e36 


0.0046 


0.13 


2.41 


1.0 


0.75 




0.57 


0.012 


5.5e36 


0.0065 


0.20 


2.41 


1.4 


1.1 


1 


0.40 


















0.57 


0.0016 


2.7e37 


0.032 


2.5 


2.33 


0.63 


1.7 


oo 


0.20 


0.0078 


6.8e37 


0.080 


5.5 


2.28 


0.52 


1.9 


10 


0.20 


0.0090 


1.7e37 


0.020 


1.3 


2.35 


0.50 


1.2 




0.28 


0.0014 


2.4e37 


0.029 


2.1 


2.35 


0.66 


1.7 


1 


0.40 


9.7E-4 


4.6e36 


0.0054 


0.35 


2.38 


0.72 


0.64 


oo 


0.20 


0.021 


9.3e36 


0.011 


0.59 


2.32 


0.52 


0.55 




0.28 


0.0061 


1.3e37 


0.016 


0.92 


2.32 


0.69 


0.78 


10 


0.28 


0.0083 


3.5e36 


0.0041 


0.23 


2.40 


0.67 


0.52 



j synch 
^2-10 
ergs s~' 



= 10 9 

ynch a Synch b 



,.syncl 
J2-10 



synch 
10 



/GeV° /VHE 



50 



2.0 



2.2 



2.4 



100 



2.0 



2.2 



150 



2.4 



2.0 
2.2 



2.4 



0.039 
0.025 
0.040 
0.027 e ' f 
0.015 
0.068 
0.035 
0.0064 
0.058 
0.039 
0.012 
0.011 
0.0024 
0.0085 
0.0094 
0.0016 
0.0012 
0.020 
0.0061 
0.0084 



0.18 


3.6e38 


0.43 


3.9 


2.08 


1.0 


2.0 


0.18 


3.6e38 


0.43 


3.9 


2.08 


1.0 


2.0 


0.17 


3.5e38 


0.41 


3.8 


2.03 


1.1 


2.0 


0.14 


2.9e38 


0.34 


3.5 


2.03 


1.3 


2.0 


0.43 


5.1e37 


0.060 


0.33 


2.27 


0.95 


0.51 


0.41 


5.0e37 


0.059 


0.34 


2.27 


1.0 


0.61 


0.34 e 


4.5e37 


0.053 


0.35 


2.27 


1.2 


0.96 


0.23 


4.1e37 


0.048 


0.38 


2.27 


1.9 


1.9 


0.44 


5.1e37 


0.061 


0.33 


2.20 


0.95 


0.50 


0.42 


4.9e37 


0.058 


0.33 


2.20 


1.0 


0.56 


0.31 e 


3.8e37 


0.045 


0.31 


2.20 


1.3 


0.88 


0.24 


3.1e37 


0.037 


0.29 


2.20 


1.9 


1.3 


0.45 


5.3e37 


0.062 


0.33 


2.12 


0.95 


0.50 


0.44 


5.1e37 


0.061 


0.33 


2.12 


1.0 


0.51 


0.27 


3.2e37 


0.038 


0.28 


2.12 


1.6 


0.61 


0.048 


1.0e39 


1.2 


49 


2.00 


0.51 


2.0 


0.038 


1.4e38 


0.16 


5.4 


2.19 


0.54 


1.6 



l.le38 
8.5e37 
l.le38 
7.2e37 
4.3e37 
1.9e37 
2.1e37 
2.3e37 
1.2e37 
1.0e37 
7.4e36 
2.6e38 
8.1e37 
9.2e37 
4.4e37 
2.9e37 
8.0e36 
1.3e37 
1.4e37 
5.0e36 



0.13 
0.10 
0.12 
0.085 
0.051 
0.022 
0.025 
0.027 
0.014 
0.012 
0.0087 
0.31 
0.095 
0.11 
0.051 
0.034 
0.0094 
0.015 
0.017 
0.0059 



4.3 
3.5 
3.9 
2.8 
1.7 
0.56 
0.70 
0.84 
0.37 
0.33 
0.26 
23 
7.4 
7.4 
3.5 
2.4 
0.61 
0.83 
1.0 
0.34 



2.14 
2.14 
2.08 
2.08 
2.08 
2.29 
2.29 
2.29 
2.26 
2.26 
2.26 
2.01 
2.01 
2.20 
2.15 
2.15 
2.10 
2.28 
2.28 
2.27 



0.53 
0.80 
0.56 
0.85 
1.1 
0.62 
1.0 
1.4 
0.78 
1.0 
1.4 
0.50 
0.63 
0.52 
0.50 
0.66 
0.72 
0.52 
0.69 
0.67 



1.1 
1.9 

0.59 
0.82 
1.0 
0.61 
1.2 
1.6 
0.57 
0.78 
1.1 
1.4 
1.7 
1.9 
1.2 
1.7 
0.64 
0.56 
0.78 
0.53 



a Fraction of the observed diffuse 2-10 keV X-ray luminosity (from Bauer et al. 2008) from that synchrotron emission accounts for. 

Ratio of the 2-10 keV synchrotron and IC luminosities. 
c Ratio of the predicted 0.3 - 10 GeV 7-ray emission with the flux detected by F<?rm/-LAT in Abdo et al. 1 2010b). 

Ratio of the predicted >220 GeV 7-ray emission with the observed TeV flux observed from HESS in Acero et al. 1 2009). 

e This model satisfies our constraints (§ !3.3t . While 7? is not the minimum or maximum allowed or the closest to reproducing the GeV luminosity for = 10 9 , this is true for 

7^"™ = 10 6 so we list it anyway, 
f 

Chosen as our fiducial model. 



TABLE 4 
M82 X-RAY LUMINOSITIES 



B p r diff (3GeV) r, = 10° 7=' = 10 y 

j. j synch r synch a Synch h x^synch /■ c r d f T synch ,-synch a Synch u y^synch r c J- d 

? L 2-& h-t, IC - 1 2-8 JGeV /VHE C ^ 2 -8 /2-8 Tc - 1 2-8 JOcV /VHE 

fiG (Myr) ergs s ergs s~' 

50 2.0 00 0.013 0.067 5.4e37 0.012 0.29 2.25 1.8 0.63 

0.035 0.066 1.5e38 0.035 0.83 2.25 1.9 1.8 
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TAB LE 4 — Continued 



B 


p fdfflO GeV) 


'/ 






prim _ 1q6 
/max 1V/ 






Pnm _ 1Q 9 
/max lu 










. synch 
^2-8 


rsynch a Synch b -psynch r c r d 
■'2-8 ~lC~ i 2-8 /GeV / VHE 




. synch 
^2-8 


^synch a Synch b pSynch /■ c r d 
^2-8 "Tc~ 1 2-8 /GeV /VHE 




(Myr) 






ergs S 






ergs S 











0.035 


0.066 


2.5e37 


0.0056 


0.14 


2.31 








0.071 


0.065 


5.0e37 


0.011 


0.27 


2.31 




2.2 


oo 


0.0088 
0.025 


















0.035 


0.10 


3.3e37 


0.0075 


0.13 


2.28 








0.10 


0.094 


9.4e37 


0.021 


0.36 


2.28 






10 


0.0088 
0.071 


















0.10 


0.096 


1.6e37 


0.0036 


0.061 


2.34 








0.14 


0.091 


2.2e37 


0.0051 


0.087 


2.34 






1 


0.0088 
0.14 














2.4 


oo 


0.050 


















0.14 


0.22 


1.7e37 


0.0040 


0.047 


2.31 






10 


0.14 












100 


2.0 


oo 


0.0088 


















0.013 


0.023 


8.3e37 


0.019 


1.3 


2.24 








0.035 


0.022 


2.4e38 


0.054 


3.7 


2.24 






10 


0.0088 


















0.035 


0.022 


4.0e37 


0.0092 


0.63 


2.30 








0.10 


0.018 


l.le38 


0.026 


1.8 


2.30 






1 


0.0088 


















0.20 


0.018 


2.7e37 


0.0062 


0.41 


2.32 








0.57 


0.0047 


7.8e37 


0.018 


1.2 


2.32 




2.2 


oo 


0.025 


















0.035 


0.031 


5.3e37 


0.012 


0.63 


2.27 








0.10 


0.022 


1.5e38 


0.034 


1.9 


2.27 






10 


0.050 


















0.10 


0.023 


2.6e37 


0.0060 


0.32 


2.34 








0.28 


7.5E-4 


7.5e37 


0.017 


1.1 


2.34 






1 


0.40 














2.4 


oo 


0.14 


0.042 


2.9e37 


0.0065 


0.28 


2.31 








0.28 


0.0081 


5.7e37 


0.013 


0.64 


2.31 


150 


2.0 


10 
1 


0.10 
0.14 


0.0075 


1.4e38 


0.032 


4.1 


2.28 








0.20 


0.0067 


3.4e37 


0.0077 


0.93 


2.30 








0.28 


0.0036 


4.8e37 


0.011 


1.3 


2.30 




2.2 


oo 


0.071 


0.0099 


1.3e38 


0.029 


3.0 


2.25 








0.10 


0.0061 


1.8e38 


0.041 


4.5 


2.25 






10 


0.10 


0.0071 


3.3e37 


0.0075 


0.80 


2.33 








0.14 


0.0020 


4.6e37 


0.011 


1.2 


2.33 




2.4 


oo 


0.14 


0.0060 


3.5e37 


0.0079 


0.74 


2.29 


200 


2.0 


1 


0.20 


0.0021 


3.8e37 


0.0087 


1.6 


2.28 




2.2 


10 


0.10 


7.2E-4 


3.7e37 


0.0085 


1.5 


2.31 



1.8 0.61 
2.0 1.2 







0.12 


7.1e38 


0.16 


2.5 


2.17 


1.4 


1.8 






0.11 


6.9e38 


0.16 


2.5 


2.17 


1.5 


2.0 


1.4 


0.61 
















1.8 


1.7 




















0.12 


7.1e38 


0.16 


2.5 


2.07 


1.4 


1.7 






0.11 


6.3e38 


0.14 


2.4 


2.07 


1.7 


1.9 


1.8 


0.68 
















2.0 


0.96 




















0.13 


7.5e38 


0.17 


2.6 


2.01 


1.5 


1.7 






0.11 


6.3e38 


0.14 


2.3 


2.01 


1.8 


1.6 






0.28 


1.0e38 


0.023 


0.22 


2.26 


1.5 


0.52 


1.9 


0.56 


0.21 


9.0e37 


0.020 


0.24 


2.26 


1.9 


0.81 






0.22 


7.9e37 


0.018 


0.20 


2.19 


1.9 


0.50 






0.037 


2.3e39 


0.52 


35 


2.12 


0.62 


1.9 


0.59 


0.55 
















0.67 


1.6 




















0.037 


2.2e39 


0.51 


35 


2.03 


0.61 


1.7 


0.65 


0.57 


0.034 


2.1e39 


0.48 


33 


2.03 


0.69 


2.0 


0.84 


1.6 




















0.039 


2.4e39 


0.54 


36 


1.98 


0.61 


1.6 


0.93 


0.54 


0.028 


1.7e39 


0.39 


26 


1.98 


0.96 


1.7 


1.6 


1.5 


0.0071 


5.1e38 


0.12 


7.6 


1.98 


1.6 


1.8 






0.034 


2.9e38 


0.067 


3.5 


2.20 


0.52 


0.65 


0.57 


0.56 


0.033 e 


3.0e38 


0.068 


3.6 


2.20 


0.58 


0.80 


0.94 


1.6 


0.024 


3.3e38 


0.074 


4.1 


2.20 


0.94 


1.8 






0.032 


2.5e38 


0.057 


3.0 


2.13 


0.64 


0.55 


0.89 


0.65 


0.025 


2.1e38 


0.049 


2.6 


2.13 


0.90 


0.83 


1.8 


1.8 


9.4E-4 


8.2e37 


0.019 


1.2 


2.13 


1.8 


1.8 






0.0064 


6.1e37 


0.014 


0.73 


2.07 


1.8 


0.50 


1.1 


0.52 


0.042 


4.8e37 


0.011 


0.47 


2.28 


1.1 


0.55 


2.0 


1.0 


0.0081 


6.1e37 


0.014 


0.68 


2.28 


2.0 


1.1 


0.56 


1.5 


0.012 


9.0e38 


0.20 


26 


2.05 


0.57 


1.8 






0.014 


9.2e38 


0.21 


25 


2.00 


0.55 


0.70 


0.65 


0.52 


0.010 e 


7.0e38 


0.16 


19 


2.00 


0.65 


0.77 


0.81 


0.74 


0.0056 


4.1e38 


0.093 


11 


2.00 


0.81 


0.87 


0.58 


1.1 


0.011 


2.1e38 


0.049 


5.1 


2.20 


0.58 


1.1 


0.74 


1.5 


0.0064 


2.3e38 


0.053 


5.9 


2.20 


0.74 


1.5 


0.70 


0.63 


0.0077 f 


9.7e37 


0.022 


2.3 


2.15 


0.70 


0.66 


0.91 


0.89 


0.0023 


6.5e37 


0.015 


1.7 


2.15 


0.91 


0.90 


0.97 


0.51 


0.0061 


3.8e37 


0.0087 


0.81 


2.27 


0.97 


0.51 


0.54 


0.52 


0.0032 


2.5e38 


0.057 


10 


2.01 


0.54 


0.56 


0.62 


0.62 


7.0E-4 


4.4e37 


0.0099 


1.7 


2.16 


0.62 


0.62 



3 Fraction of the observed diffuse 2-8 keV X-ray luminosity (from Strickland & Heckman 2007) from that synchrotron emission accounts for. 

Ratio of the 2-8 keV synchrotron and IC luminosities. 
C Ratio of the predicted 0.3 - 10 GeV 7-ray emission with the flux detected by Fermi-LAT in Abdo et al. 1 2010b). 
^ Ratio of the predicted >700 GeV 7-ray emission with the observed TeV flux observed from HESS in Acciari et al. 1 2009). 

6 This model satisfies our constraints f§ |3.3| . While r/ is not the minimum or maximum allowed or the closest to reproducing the GeV luminosity for = 10 9 , this is true for 

7max = so we list it anyway, 
f 

Chosen as our fiducial model. 

TABLE 5 
ARP 220 X-RAY LUMINOSITIES 



B p , dlff (3GeV) 7n P , r r = 10 6 7 mT = 10 9 

j- .synch ,.synch a Synch b psynch r c r d a .synch ,.synch a Synch b psynch r c 

S L 2-I0 J2-10 IC 2-10 J GcV /VHE ? ^2-10 J2-10 IC 2-10 J GeV 

/tG (Myr) ergs s -1 ergs s 



West nucleus 



500 2.0 


oo 


0.81 


2.6e39 


0.065 


0.0416 


1.90 


0.60 


0.016 


1.3 


1.2e41 


2.9 


1.9 


1.85 


0.67 


0.23 




10 


0.81 


2.4e39 


0.061 


0.0388 


1.90 


0.60 


0.016 


1.3 


1.2e41 


2.9 


1.9 


1.85 


0.67 


0.23 




1 


0.81 


1.5e39 


0.039 


0.0246 


1.93 


0.60 


0.012 


1.3 


1.2e41 


2.9 


1.9 


1.84 


0.66 


0.23 


2.2 


oo 


1.3 


6.0e38 


0.015 


0.0072 


1.94 


0.49 


0.0060 


1.3 


1.7e40 


0.41 


0.20 


1.91 


0.50 


0.044 




10 


1.3 


5.6e38 


0.014 


0.0068 


1.94 


0.49 


0.0058 


1.3 


1.7e40 


0.41 


0.20 


1.91 


0.50 


0.044 




1 


1.3 


3.6e38 


0.0091 


0.0044 


1.97 


0.49 


0.0047 


1.3 


1.6e40 


0.41 


0.20 


1.92 


0.50 


0.043 


2.4 


oo 


2.8 


8.5e37 


0.0021 


8E-4 


1.98 


0.46 


0.0015 


2.8 


2.3e39 


0.058 


0.021 


1.97 


0.46 


0.0085 




10 


2.8 


8.0e37 


0.0020 


7E-4 


1.99 


0.46 


0.0014 


2.8 


2.3e39 


0.058 


0.021 


1.98 


0.46 


0.0085 




1 


2.8 


5.3e37 


0.0013 


5E-4 


2.02 


0.46 


0.0012 


2.8 


2.3e39 


0.058 


0.021 


2.00 


0.46 


0.0082 
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TABLE 5 — Continued 

B p tof (3GeV) 7 !T= 10 6 iST = 10 9 











-^2-10 


/■synch y 
•/2-10 


Synch p 


T^synch 
* 2-10 


/GeV C 


/vHE d 




synch 
^2-10 


f*y^"A 
J2-W 


Synch ^ 


T-i.synch 
1 2-10 


/GeV C 


/vHE d 






(Myr) 




ergs s -1 














ergs s _1 












1000 


2.0 


oo 


0.25 


6.5e39 


0.16 


0.33 


1.97 


0.19 


0.015 


0.38 


7.3e40 


1.8 


3.7 


1.93 


0.21 


0.061 






10 


0.25 


6.1e39 


0.15 


0.31 


1.97 


0.19 


0.014 


0.38 


7.3e40 


1.8 


3.7 


1.93 


0.21 


0.060 






1 


0.25 


3.9e39 


0.098 


0.20 


2.00 


0.19 


0.011 


0.38 


7.1e40 


1.8 


3.6 


1.93 


0.21 


0.057 




2.2 


oo 


0.35 


1.6e39 


0.039 


0.063 


2.01 


0.16 


0.0056 


0.38 


1.1 e40 


0.28 


0.45 


1.99 


0.16 


0.014 






10 


0.35 


1.5e39 


0.037 


0.060 


2.01 


0.16 


0.0054 


0.38 


l.le40 


0.27 


0.45 


1.99 


0.16 


0.014 






1 


0.36 


9.7e38 


0.024 


0.040 


2.04 


0.16 


0.0045 


0.38 


l.le40 


0.26 


0.43 


2.01 


0.16 


0.013 




2.4 


oo 


0.75 


2.3e38 


0.0058 


0.0075 


2.06 


0.15 


0.0014 


0.76 


1.6e39 


0.041 


0.053 


2.05 


0.15 


0.0029 






10 


0.75 


2.2e38 


0.0054 


0.0071 


2.07 


0.15 


0.0014 


0.76 


1.6e39 


0.041 


0.052 


2.06 


0.15 


0.0028 






1 


0.76 


1.5e38 


0.0037 


0.0048 


2.10 


0.15 


0.0012 


0.76 


1.6e39 


0.039 


0.050 


2.08 


0.15 


0.0026 


2000 


2.0 


oo 


0.077 


1.3e40 


0.31 


1.9 


2.04 


0.069 


0.012 


0.12 


4.6e40 


1.1 


6.9 


2.01 


0.073 


0.020 






10 


0.077 


1.2e40 


0.29 


1.8 


2.04 


0.069 


0.012 


0.12 


4.5e40 


1.1 


6.8 


2.01 


0.073 


0.020 






1 


0.077 


7.9e39 


0.20 


1.2 


2.07 


0.068 


0.0096 


0.12 


4.1e40 


1.0 


6.2 


2.02 


0.048 


0.00/8 




O O 

1.1 


oo 


0.097 


3.2e39 


0.079 


n a i 
0.41 


o nn 
2.09 


0.064 


n nn/in 
U.U049 


0. 10 


7.6e39 


0.19 


1.0 


o no 
2.0/ 


0.064 


0.0062 






1 n 
10 


0.097 


3.0e39 


0.075 


0.39 


o nn 
2.09 


0.064 


n r\t\A q 


0. 10 


7.5e39 


0. 19 


0.98 


o no 
2.Uo 


0.065 


0.0061 






1 


0.09 / 


2.0e39 


0.051 


0.26 


2.12 


0.063 


n dd/t n 
U.U040 


0. 10 


6.6e39 


0.17 


0.86 


2. 10 


0.064 


0.0053 




/I 
Z.4 


oo 


A on 
U.ZU 


4.ye3o 


\J\)YL 


U.UD4 


IK 
Z. 1 J 


n nAi 
U.UD1 


n nm o 
U.UU13 


n on 
U.ZU 


1 OaOO 

i.ze3y 


n nm 
U.U3U 


n 1 1 
U. 13 


O 1 A 

Z.14 


n nA i 
U.Uol 


n nn i 5 
U.UU1 j 






1 n 


n on 
U.ZU 


A OaOC 

4. /e3o 


U.Ulz 


n n^o 
U.UjZ 


1 5 
Z. 1 J 


n no.< 
U.U3J 


n fiftl o 
U.UU13 


n on 
U.zU 


1 OaOO 

i.ze3y 


n non 

u.uzy 


nil 
U. 13 


O 1 a 
Z.l J 


n n£ i 
U.Uol 


n nn i 5 
U.UU1 j 






i 


n on 
U.ZU 


3.3e3o 


U.UUoZ 


n noA 
U.U30 


T in 

i. iy 


n n£i 
U.UD 1 


n nm i 
U.UU1 1 


n on 
U.ZU 


i.ue3y 


n 

U.UZO 


n 1 1 
U. 1 1 


O 1 o 

z.l / 


n n/s i 
U.Uol 


n nn i o 
U.UU13 


A (1(1(1 

4000 


O A 

2.0 


oo 


n noo 
0.023 


1.8e40 


A AC 


O 
/.5 


o no 
2.0/ 


n non 
0.030 


n m n 
0.010 


0.03O 


o 1 a/in 
3.1e40 


n on 

u. /y 


1 a 
13 


O (15 

2.05 


n non 
0.030 


n (i i 1 
U.U1 1 






1 n 
10 


n noo 
0.023 


1 Oa/1(l 

1 . / e40 


A A A 

U.44 


o a 
/.3 


o no 
2.0/ 


n non 
0.030 


n nnnn 
0.0099 


n no£ 
0.03O 


o n^/in 
3.0e40 


n 76. 
U. /D 


13 


O (15 

z.Uj 


n non 
0.030 


n m 1 
0.01 1 






1 


0.024 


1.2e40 


0.30 


5.0 


2.10 


0.029 


0.0082 


0.036 


2.5e40 


0.63 


11 


2.07 


0.030 


0.0091 




2.2 


oo 


0.020 


4.9e39 


0.12 


1.9 


2.13 


0.036 


0.0043 


0.021 


6.2e39 


0.16 


2.4 


2.12 


0.036 


0.0044 






10 


0.020 


4.6e39 


0.12 


1.8 


2.13 


0.036 


0.0042 


0.021 f 


6.0e39 


0.15 


2.3 


2.12 


0.036 


0.0043 






1 


0.020 


3.3e39 


0.082 


1.3 


2.16 


0.035 


0.0036 


0.022 


4.6e39 


0.12 


1.8 


2.15 


0.035 


0.0037 




2.4 


oo 


0.039 


8.0e38 


0.020 


0.28 


2.19 


0.037 


0.0012 


0.039 


1.0e39 


0.026 


0.36 


2.19 


0.037 


0.0012 






10 


0.039 


8.0e38 


0.020 


0.28 


2.19 


0.037 


0.0012 


0.040 


9.8e38 


0.025 


0.35 


2.19 


0.037 


0.0012 






1 


0.040 


5.5e38 


0.014 


0.20 


2.23 


0.037 


0.0010 


0.040 


7.7e38 


0.019 


0.27 


2.22 


0.037 


0.0010 


8000 


2.0 


oo 


0.0057 


2.2e40 


0.56 


26 


2.06 


0.018 


0.0092 


0.0088 


2.6e40 


0.65 


30 


2.05 


0.018 


0.0093 






10 


0.0057 


2.1e40 


0.53 


24 


2.06 


0.018 


0.0090 


0.0088 


2.5e40 


0.62 


28 


2.05 


0.018 


0.0091 






1 


0.0059 


1.5e40 


0.37 


17 


2.09 


0.018 


0.0075 


0.0092 


1.9e40 


0.46 


21 


2.07 


0.018 


0.0076 


16000 


2.0 


1 


1.8E-4 


1.6e40 


0.41 


57 


2.06 


0.015 


0.0073 


2.7E-4 


1.7e40 


0.41 


57 


2.04 


0.015 


0.0073 


East nucleus 


500 


2.0 


oo 


0.61 


2.6e39 


0.17 


0.055 


1.90 


0.46 


0.016 


0.95 


8.9e40 


5.9 


1.9 


1.85 


0.51 


0.18 






10 


0.61 


2.4e39 


0.16 


0.051 


1.90 


0.46 


0.016 


0.95 


8.9e40 


5.9 


1.9 


1.85 


0.51 


0.18 






1 


0.61 


1.5e39 


0.10 


0.032 


1.93 


0.46 


0.012 


0.95 


8.8e40 


5.9 


1.9 


1.84 


0.51 


0.17 




2.2 


oo 


0.95 


6.0e38 


0.040 


0.0095 


1.94 


0.38 


0.0060 


1.0 


1.3e40 


0.84 


0.20 


1.91 


0.39 


0.035 






10 


0.95 


5.6e38 


0.038 


0.0089 


1.94 


0.38 


0.0058 


1.0 


1.3e40 


0.84 


0.20 


1.91 


0.39 


0.035 






1 


0.95 


3.6e38 


0.024 


0.0057 


1.97 


0.37 


0.0047 


1.0 


1.2e40 


0.83 


0.20 


1.92 


0.39 


0.034 




2.4 


oo 


2.1 


8.5e37 


0.0057 


0.0010 


1.98 


0.36 


0.0015 


2.1 


1.8e39 


0.12 


0.021 


1.97 


0.36 


0.0068 






10 


2.1 


8.0e37 


0.0054 


9E-4 


1.99 


0.36 


0.0014 


2.1 


1.8e39 


0.12 


0.021 


1.98 


0.36 


0.0067 






1 


2.1 


5.3e37 


0.0035 


6E-4 


2.02 


0.36 


0.0012 


2.1 


1.8e39 


0.12 


0.021 


2.00 


0.36 


0.0065 


1000 


2.0 


oo 


0.18 


6.5e39 


0.43 


0.43 


1.97 


0.15 


0.015 


0.29 


5.6e40 


3.7 


3.8 


1.93 


0.16 


0.049 






10 


0.18 


6.1e39 


0.40 


0.41 


1.97 


0.15 


0.014 


0.29 


5.6e40 


3.7 


3.8 


1.93 


0.16 


0.049 






1 


0.18 


3.9e39 


0.26 


0.26 


2.00 


0.15 


0.011 


0.29 


5.4e40 


3.6 


3.6 


1.93 


0.16 


0.046 




2.2 


oo 


0.26 


1.6e39 


0.10 


0.084 


2.01 


0.13 


0.0056 


0.28 


8.6e39 


0.57 


0.46 


1.99 


0.13 


0.012 






10 


0.26 


1.5e39 


0.098 


0.079 


2.01 


0.13 


0.0054 


0.28 


8.5e39 


0.57 


0.46 


1.99 


0.13 


0.011 






1 


0.26 


9.7e38 


0.065 


0.052 


2.04 


0.13 


0.0045 


0.28 


8.0e39 


0.53 


0.43 


2.01 


0.13 


0.010 




2.4 


oo 


0.56 


2.3e38 


0.015 


0.0099 


2.06 


0.12 


0.0014 


0.56 


1.3e39 


0.085 


0.054 


2.05 


0.12 


0.0025 






10 


0.56 


2.2e38 


0.015 


0.0093 


2.07 


0.12 


0.0014 


0.56 


1.3e39 


0.084 


0.054 


2.06 


0.12 


0.0025 






1 


0.56 


1.5e38 


0.0099 


0.0063 


2.10 


0.12 


0.0012 


0.56 


1.2e39 


0.079 


0.051 


2.08 


0.12 


0.0023 


2000 


2.0 


oo 


0.054 


1.3e40 


0.83 


2.5 


2.04 


0.056 


0.012 


0.084 


3.6e40 


2.4 


7.2 


2.01 


0.059 


0.018 






10 


0.054 


1.2e40 


0.79 


2.4 


2.04 


0.056 


0.012 


0.084 


3.5e40 


2.3 


7.1 


2.01 


0.059 


0.017 






1 


0.055 


7.9e39 


0.52 


1.6 


2.07 


0.055 


0.0096 


0.085 


3.1e40 


2.1 


6.3 


2.02 


0.058 


0.015 




2.2 


oo 


0.066 


3.2e39 


0.21 


0.54 


2.09 


0.054 


0.0049 


0.070 


6.3e39 


0.42 


1.1 


2.07 


0.055 


0.0058 






10 


0.066 


3.0e39 


0.20 


0.51 


2.09 


0.054 


0.0048 


0.070 


6.1e39 


0.41 


1.0 


2.08 


0.055 


0.0057 






1 


0.066 


2.0e39 


0.14 


0.35 


2.12 


0.054 


0.0040 


0.071 


5.2e39 


0.34 


0.89 


2.10 


0.054 


0.0049 




2.4 


oo 


0.14 


4.9e38 


0.033 


0.072 


2.15 


0.053 


0.0013 


0.14 


9.8e38 


0.065 


0.14 


2.14 


0.053 


0.0015 






10 


0.14 


4.7e38 


0.031 


0.069 


2.15 


0.053 


0.0013 


0.14 


9.5e38 


0.064 


0.14 


2.15 


0.053 


0.0014 






1 


0.14 


3.3e38 


0.022 


0.048 


2.19 


0.053 


0.0011 


0.14 


8.2e38 


0.054 


0.12 


2.17 


0.053 


0.0013 


4000 


2.0 


oo 


0.014 


1.8e40 


1.2 


10 


2.07 


0.026 


0.010 


0.021 


2.6e40 


1.7 


15 


2.05 


0.026 


0.011 






10 


0.014 


1.7e40 


1.2 


9.7 


2.07 


0.026 


0.0099 


0.021 


2.5e40 


1.7 


14 


2.05 


0.026 


0.011 






1 


0.014 


1.2e40 


0.80 


6.7 


2.10 


0.025 


0.0082 


0.022 


2.0e40 


1.3 


11 


2.07 


0.026 


0.0087 




2.2 


oo 


0.0077 


4.9e39 


0.33 


2.5 


2.13 


0.033 


0.0043 


0.0082 


5.4e39 


0.36 


2.8 


2.12 


0.033 


0.0044 






10 


0.0078 


4.6e39 


0.31 


2.4 


2.13 


0.033 


0.0042 


().()083 f 


5.2e39 


0.34 


2.7 


2.12 


0.033 


0.0043 






1 


0.0082 


3.3e39 


0.22 


1.7 


2.16 


0.033 


0.0036 


0.0087 


3.8e39 


0.25 


2.0 


2.15 


0.033 


0.0037 




2.4 


oo 


0.017 


8.0e38 


0.054 


0.38 


2.19 


0.035 


0.0012 


0.017 


9.0e38 


0.060 


0.43 


2.19 


0.035 


0.0012 






10 


0.017 


7.7e38 


0.051 


0.37 


2.20 


0.035 


0.0012 


0.017 


8.6e38 


0.057 


0.41 


2.19 


0.035 


0.0012 






1 


0.017 


5.5e38 


0.037 


0.26 


2.23 


0.034 


0.0010 


0.018 


6.5e38 


0.043 


0.31 


2.22 


0.034 


0.0010 


8000 


2.0 


oo 


0.00054 


2.2e40 


1.5 


34 


2.06 


0.017 


0.0092 


8.4E-4 


2.3e40 


1.5 


35 


2.05 


0.017 


0.0092 






10 


0.00057 


2.1e40 


1.4 


32 


2.06 


0.017 


0.0090 


8.9E-4 


2.2e40 


1.4 


33 


2.05 


0.017 


0.0090 






1 


0.00083 


1.5e40 


0.99 


23 


2.09 


0.017 


0.0075 


0.0013 


1.5e40 


1.0 


23 


2.07 


0.017 


0.0076 
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TABLE 5 — Continued 



B 


p t m Q GeV) 






prim _ 1q6 
/max 1 u 






prim _ 1Q 9 
/max 1 u 






e 


. synch 
^2-10 


r synch a Synch h -i-iSynch /■ c r d 
J2-1Q ~lC~ 1 2-10 /GeV /VHE 




j synch 
^2-10 


f synch a Synch b -psynch /■ c x d 
/2-10 ~lC~ 1 2-10 /GeV /VHE 


f iG 


(Myr) 




ergs s -1 






ergs s 





a Fraction of the observed diffuse 2-10 keV X-ray luminosity (from Clements et al. 2002) from that synchrotron emission accounts for. We assume that Arp 220 X-l is the west 

nucleus, and Arp 220 X-4 is the east nucleus. 

k Ratio of the 2 - 10 keV synchrotron and IC luminosities. 

C Ratio of the predicted > 100 MeV 7-ray emission with the Fermi-hAH one-year catalog source sensitivity at high Galactic latitude 1 Abdo et al. 20 1 0b ). 
^ Ratio of the predicted 0.36 - 1 .8 TeV 7-ray emission with the upper limits on the flux obtained with MAGIC in Alb ert et alj |2007f) . 
6 Chosen as our fiducial model. 



